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abstract. We prove existence of a special class of solutions to the (elliptic) Nonlinear Schrodinger 
Equation — e 2 At/j + V(x)ip — \ip\ p ~ 1 %p, on a manifold or in the Euclidean space. Here V represents the 
potential, p an exponent greater than 1 and e a small parameter corresponding to the Planck constant. As 
e tends to zero (namely in the scmiclassical limit) we prove existence of complex- valued solutions which 
concentrate along closed curves, and whose phase is highly oscillatory. Physically, these solutions carry 
quantum-mechanical momentum along the limit curves. In this first part we provide the characterization 
of the limit set, with natural stationarity and non-degeneracy conditions. We then construct an approx- 
imate solution up to order e 2 , showing that these conditions appear naturally in a Taylor expansion of 
the equation in powers of e. Based on these, an existence result will be proved in the second part |i0] . 
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1 Introduction 

In this paper we are concerned with concentration phenomena for solutions of the singularly-perturbed 
elliptic equation 

(NLS £ ) -e 2 A g tP + V(x)4;= \ip\ p - l ip on M, 

where M is an n-dimensional compact manifold (or the flat Euclidean space R"), V a smooth positive 
function on M satisfying the properties 

(1) 0<V 1 <V<V 2 ; \\V\\ C 3<V 3 , 

(for some fixed constants Vi, 1 ^,!^) tp a complex- valued function, e > a small parameter and p an 
exponent greater than 1. Here A g stands for the Laplace-Beltrami operator on (M,g). 
( NLS S \ arises from the study of the Nonlinear Schrodinger Equation 

(2) ih^- = -h 2 A4> + V(x)<ip - \4>\ p ~ 1 ^ on M x [0, +oo), 

at 
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where ip — ip(x,t) is the wave function, V(x) a potential, and h the Planck constant. A special class 
of solutions to ([2| is constituted by the functions whose dependence on the variables x and t is of the 
form ij)(x,t) — e~ 1 ^ ^p(x). Such solutions are called standing waves and up to substituting V(x) with 
V(x) — lu, they give rise to solutions of {NLS E ), for e = h. 

An interesting case is the semiclassical limit e — > 0, where one should expect to recover the Newton 
law of classical mechanics. In particular, near stationary points of the potential, one is lead to search 
highly concentrated solutions, which could mimic point-particles at rest. 

In recent years, a lot of attention has been devoted to the study of the above problem: one of the 
first results in this direction is due to Floer and Weinstein in [2B], where the case of M — K and p = 3 
is considered, and where existence of solutions highly concentrated near critical points of V has been 
proved. This result has then been extended by Oh, [5T], to the case of R n for arbitrary n, provided 
1 < P < The profile of these solutions is given by the ground state U Xo (namely the solution with 

minimal energy, which is real-valued, everywhere positive and can be assumed radial) of the following 
limit equation 



(3) 



Au + V(x Q )u = u p 



where xq is the concentration point. The solutions u e obtained in the aforementioned papers behave 
qualitatively like u s (x) ~ U Xg ( x ~ x ° \ as e tends to zero, and since U Xo decays exponentially to zero at 
infinity, u e vanishes rapidly away from x$. 

Two comments are in order: first of all the criticality of V at xq is a necessary condition for such a 
behavior, as shown in |55j . Secondly, as pointed out in |17j . also the upper bound p < is required 
for having solutions concentrating at points: indeed the well-known Pohozaev's identity imposes this 
restriction for having existence of solutions to ^ tending to zero at infinity. 

The above existence results have been extended in several directions, including the construction of 
solutions with multiple peaks, the case of degenerate potentials, potentials tending to zero at infinity and 
more general nonlincarities. We refer the interested reader for example to the (incomplete) list of works 
DP, 0, E], [6], 0, Ej, [32], [22], [28], [33] and to the bibliog raphies therein. 

We also mention the mathematical similarities between ( NLS £ I and the following problem 



—e 2 Au + u = u p 



(Pe) 



du 



= 



du 

u > 



in f2, 
on d£l, 
in f2, 



where ft is a smooth bounded domain of K , p > 1, and where v denotes the exterior unit normal vector 



to dil. Problem (P € \ arises in the study of some biological models, see for example |47j and references 
therein, and as (NLS e I it exhibits concentration of solutions at some points of il. Since the last equation 



is homogeneous, the location of the concentration points is determined by the geometry of the domain: if 
it occurs at the boundary, these are critical points of the mean curvature while if it occurs at the interior 
these points are (roughly) singular points for the distance function from the boundary. About this topic, 
we refer the reader to [19], ED, |29], ED, EH, [32], [35], [36], [37], 08], 09], [50], [56]. 



More recently, new types of solutions to (NLS E \ have been found, since when e tends to zero they 



do not concentrate at points, but instead at sets of higher dimension. Before stating our main result, it 
is convenient to recall the progress on this topic and to illustrate the new phenomena involved. Some 
first results in this direction were given in [12], [14] in the case of radial symmetry, and later improved in 
[1] (see also [5] for the problem in bounded domains), where necessary and sufficient conditions for the 
location of the concentration set have been given. Differently from the previous case, the limit set is not 
anymore stationary for the potential V: indeed, from heuristic considerations, the energy of a solution 
concentrated near a sphere of radius r depends both on V and on its volume, which is proportional 
to er n_1 . In [4] it was shown that the candidate radii of concentration are the critical points of the 
function r™~ 1 I/(r)p TrT_5 (the power of V in this formula arises from some scaling argument, related to 
the dependence in V(xq) of the solutions to see also Section |2j. Furthermore, no upper bound on 



2 



the exponent p is required: in fact the profile of these solutions is given by the solution of Q in K 1 , and 
in one dimension there is no restriction for the existence of entire solutions. 

Based on the above energy considerations, in |4j it is also stated a conjecture concerning concentration 
on fc-dimensional manifolds, for k = 1, . . . , n— 1: it is indeed expected that, under suitable non-degeneracy 
assumptions, the limit set should satisfy the equation 

(4) 9 k v N v = vn, with 9 i = E±l-h n - k ) 

p — 1 2 

where W N stands for the normal gradient, H the curvature vector, and the profile of the solutions at a 
point xq in the limit set should be asymptotic, in the normal directions, to the ground state of 



(5) - Au+V(x a )u = u p 



i—k 



Since the Pohozaev identity implies p < "_j?^2 f° r ^ ne existence of non trivial solutions, the latter 
condition is expected to be a natural one for dealing with this phenomenon. 



Actually, concerning (P e ) another conjecture has been previously stated, asserting existence of solu- 
tions concentrating at sets of positive dimension. Concerning the latter problem, starting from the paper 
[43], there has been some progress in the general setting (without symmetry assumptions), and after the 
works [39., |42j . |44j . existence in now known for arbitrary dimension and codimcnsion. About problem 
( NLS e \ , the conjecture in [3] has been verified in [21] for n — 2 and k = 1. Some other (and related) 



results, under some reduced symmetry assumptions (as cylindrical or similar) have been given in [13] . 
[2D] . [3PJ, [S3- Specifically we mention the note [ST]: instead of it is considered there the nonlinear 
wave equation 

(6) - -4- = -A$ + V(x)<tp - M^i! on M x [0,+oo) 

ot 

for p = 3. It has been proved in [5] and [TUJ that 

d 2 ib 

has solutions which remain concentrated near elliptic closed geodesies in M for long periods of time, 
but which eventually drift away when t — > oo. A. Weinstein in |57] proved that there do exist periodic 
solutions of ^ remaining concentrated for all times, whenever p = 3, M = S 2 and V is odd. 

It is worth pointing out a major difference between the symmetric and the non-symmetric situation. 
In fact, since the ground states of ^ or ^ are of mountain-pass type (namely critical points of some 
Euler functional with Morse index equal to 1), equation ( NLS e ty becomes highly resonant. To explain the 



reason, we consider for example a real-valued function tp in R^with a radial potential. One can begin by 
finding approximate (radial) solutions of the form Ur tS (r) ~ Up (^j 1 ), where Ur is the solution of Q for 
n = 1 corresponding to V(f). Then, with a good choice of r, one can try to linearize the equation and 
find true solutions via the implicit function theorem. The linearized equation, taking ip real for simplicity, 
becomes 

-e 2 Aip + V(r)tp-pu^ e (r) p ' 1 ip in M 2 . 

Using polar coordinates (r, $) and a Fourier decomposition of ip with respect to i?, tp(r, #) = J^j el ''^' l l J j( r )> 
we see that on each component ipj acts the operator 

(7) -£ 2 V" " £ 2 l^j + V(r)^ - pur^ry- 1 ^, +^e 2 j% on [0, +oo), 

v v ' 

where L\, e (apart from the term £ 2 \^'j which is not relevant to the next discussion) represents the 



linearized equation of (NLS £ ) in one dimension near a soliton. Since one expects to deal with functions 
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which arc highly concentrated near r — r, the last term in the above formula naively increases the 
eigenvalues by a quantity of order ^£ 2 j 2 compared to those of L\ £ . 

The operator Li j£ possesses a negative eigenvalue rj £ lying between two negative constants independent 
of e (since Uy is of mountain-pass type, as explained before) and a (nearly) zero eigenvalue a £1 by the 
translation invariance of ^ in M 1 . As a consequence, the operator in ^ will possess two sequences of 
eigenvalues qualitatively of the form rjj t£ ~ i] £ + e 2 j 2 and <7j. £ ~ a £ + e 2 j 2 . This might generate two 
kinds of resonances: for small values of j, when a^ £ ~ 0, and for j of order ~, when ?7 Jj£ could be close 
to zero. A comment is in order on the corresponding eigenfunctions, which can be roughly studied with 
a separation of variables as before. The ones relative to <Tj £ (for j small) are slowly oscillating along the 
limit set, while the ones relative to the resonant rjj^'s are fast oscillating with a number of oscillations 
proportional to k ~ K 

The invcrtibility of the linearized operator will then be equivalent to having all the Uj, £ s and all the 
i]j e 's different from zero. A control on the resonant <x/, e 's can be obtained (via some careful expansions) 
from a suitable non-degeneracy condition involving the limit set and the potential V. In [4] for example, 

this can be achieved from the fact of having a non-degenerate critical point of the function r™~ 1 V(r) p^ 1 ~ 5 . 
On the other hand, the possible vanishing of some rjj tS is peculiar of this concentration behavior and more 
intrinsic, so invertibility can only be achieved by choosing suitable values of e. It is interesting to compare 



this phenomenon (which is also present in (P s )) to a result in [TS], asserting that if the Morse index of a 
family of solutions to (P £ I stays bounded as e — > 0, these must concentrate at a finite number of points. 

These formal considerations can also apply to the case of concentration near a general manifold 
(without symmetry) in higher dimension or codimension. Instead of expanding in polar coordinates one 
can use (naively) a Fourier decomposition with respect to the eigenfunctions of the Laplace-Beltrami 
operator and the normal Laplacian on the limit manifold, see |39j . If the latter has dimension k then the 
term £ 2 j 2 , by the Weyl's asymptotic formula, has to be substituted with a quantity behaving like e 2 j^ . 
We notice that in this way the average distance between two consecutive r?j, e 's (when they are close to 
zero) is of order e k so, even if we have invertibility, the distance of the spectrum to zero is (in the best 
cases) of order min{minj \r)j, £ \, min^ |oj. e |} ~ min{e 2 ,e fc }. Therefore the inverse operator is always large 
in norm. By this reason, to apply the implicit function theorem we need first to find good approximate 
solutions, with a precision depending on k, and then prove indeed that the linearized operator is invertible 
for suitable values of e. This is indeed a rather delicate issue: for reasons of brevity we do not discuss 
it here but we refer directly to |43j and |44j . Related phenomena appear in some geometric problems as 
well, dealing with the construction of surfaces with constant mean curvature, see [35], [4"5j . 

When is a radially symmetric domain and the potential V is radially symmetric, the problem is 
simpler, since working in spaces of invariant functions avoids most of the above resonances. In this case 
only finitely-many eigenvalues (depending on the dimension and the codimension of the concentration 
set) can approach zero, and localization can be determined with a finite-dimensional reduction of the 
problem. 

In this paper and in |40j we construct a new type of solutions, which concentrate along some curve 
7, and which physically carry momentum along the limit set. Differently from the solutions discussed 
before, these are complex-valued and their profile near any point xo in the image of 7 is asymptotic to 
a solution to (|3| which decays exponentially to zero away from the x n axis of W 1 and is periodic in x n . 
More precisely, we consider solutions of the form 

4>(x',x n ) = e~ lfXn U(x'), x' = (xi, . . .,x n -i), 

where / is some constant and U (x') a real function. With this choice of <f>, the function U satisfies 

(8) -AU+(f 2 + V(x )^jU=\U\ p ~ 1 U inW 1 - 1 , 

and decays to zero at infinity. Solutions to ^ can be found by considering the (real) function U satisfying 
— AU + U = U p in M n_1 (decaying to zero at infinity), and by using the scaling 

1 1 

(9) U(x')=hU(kx'), h = (/ 2 + V(x j) P " , k={f 2 + V{x )Y . 
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In the above formulas / can be taken arbitrarily, and h, k have to be chosen accordingly, depending on 
V(xo). The constant / represents the speed of the phase oscillation, and is physically related to the 
velocity of the quantum-mechanical particle represented by the wave function. 

We are aware of one result only in this direction, given in |21j , where the case of an axially-symmetric 
potential is considered, and our goal here is to treat this phenomenon in a generic situation, without 
any symmetry restriction. Some of the difficulties of such an extension were naively summarized in the 
above discussion but some new ones arise, due to the fact that the standing waves are complex- valued, 
and due to their highly oscillatory phase. In this first part we determine the concentration set and show 
that its geometric characterization appears when we construct approximate solutions to (NLS £ \, while 
a full existence result will be given in [4"0] , 

Before stating our main result we discuss how to determine the limit set: if we look for a solution i/j 



to (NLS E ) with the above profile, then it should qualitatively behave as 



(10) 



iP(s,0^e- li ^h(s)U ( ' HsK 



where s stands for the arc-length parameter of 7, and £ for a system of geodesic coordinates normal to 7. 
For having more flexibility, we chose the phase oscillation to depend on the point 7(s), while h(s),k(s) 
should satisfy 



(11) 



h(s) = ((f(s)) 2 + v( 1 (s)))^ 



m = ((f(-s)f + v( 1 (- s ))y 



which is the counterpart of ([9| for a variable potential. 



The function /(s) can be (heuristically) determined using an expansion of (NLS £ \ at order e: a 



computation performed in Subsection 2.3 (see in particular formula ( 23 1 ) shows that 
(12) f'(s) ~ Ah"(s) with a = ( " - 1 - ^ - 2, 



where A is an arbitrary constant. At this point, only the curve 7 should be determined. First of all, we 
notice that the phase should be a periodic function in the length of the curve, and therefore by ( 12 1 it is 
natural to work in the class of loops 



(13) 



r:= 7 



M periodic : A J h(s) a ds = constant j 



where s stands for the arc-length parameter on 7. Problem {NLS E \ has a variational structure, with 
Euler-Lagrange functional given by 

E e ^) = \! ( £ 2 iv s vf +n*)ivf)-^ T / 

A JM P+l J M 



For a function of the form (10), by a scaling argument (see ( 26 1 ) one has 



(14) 



h(s) e ds, 



with 9=p+l - -(p-l)(n-l), 



therefore a limit curve 7 should be a critical point of the functional J h(s) ds in the class T. With a 
direct computation, see Subsection |2.4[ one can check that the extremality condition is the following 

(15) V Ny ^ n ^P_i h p-l _ 2j[ 2 h 2a 

where, as before, W N V represents the normal gradient of V and H the curvature vector of 7. Similarly, via 
some long but straightforward calculation, one can find a natural non-degeneracy condition for stationary 
points, which is expressed by the invertibility of the operator in ( |36[ ) acting on the normal sections to 7 
(we refer the reader to Section|2]for the notation used in the formula). We notice that, since formula ( 12 1 
determines only the derivative of the phase, to obtain periodicity we need to introduce some nonlocal 



terms, see (29 1. After these preliminaries, we are in position to state our main result. 
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Theorem 1.1 Let M be a compact n- dimensional manifold, let V : M — ► R be a smooth positive function 
and let 1 < p < Let 7 be a simple closed curve in M: then there exists a positive constant Ao, 



depending on V| 7 and p for which the following holds. If < A < Ao, if"/ satisfies (151 and the operator 
in (36 1 is invertible on the normal sections 0/7, there is a sequence Sk — > such that problem (NLS £k ) 



possesses solutions Tp £k having the asymptotics in (|10[), with f satisfying (12 I 



Remark 1.2 (a) The statement of Theorem \l.l\ remains unchanged if we replace M by W l (or with an 
open manifold asymptotically Euclidean at infinity) and we assume V to be bounded between two positive 
constants and for which ||V V\\ < Ci, I = 1,2,3 for some positive constants Cj. 

(b) The restriction on the exponent p is natural, since Q admits solitons if and only if p is subcritical 
with respect to the dimension n — 1. 

(c) The smallness requirement on A is technical and we believe this condition could be relaxed. Anyway, 
for ^| < p < A should have an upper bound depending on V , to have solvability of both (|9| and 



(12 I. About this condition see Remark 2.2 and Remark 2.7 in I4(rf - 

(d) Apart from the assumption on A, Theorem \l.l\ improves the result in J21f . In fact, in addition 
to removing the symmetry condition (which is the main issue), the characterization of the limit set is 
explicit, the assumptions on V are purely local, and the upper bound on p is sharp. 



(e) The existence of solutions to (NLS £ I only for a suitable sequence — > is related to the resonance 
phenomenon described above. The result can be extended to a sequence of intervals in the parameter e 
approaching zero but, at least with our proof, we do not expect to find existence for all the epsilon's. 



Taking A — (which implies /' = 0), from (11) it follows that V = hP 1 and that condition (15) is 



equivalent to Q, so as a consequence of our result we can prove the conjecture in [3] for k 
the result in 124 . 



1, extending 



Corollary 1.3 Let M be a compact Riemannian n-dimensional manifold with metric g, let V : M — > M be 
a function satisfying (JT|) and let 1 < p < Let 7 be a simple closed curve which is a non-degenerate 

p + l _ n-l 

geodesic with respect to the weighted metric V p- 1 2 g. Then there is a sequence £k — * such that 
problem (NLS £k ) possesses real-valued solutions ip £k concentrating near 7 as j — > +00 and having the 
asymptotic behavior 

V(7(*))-: 



il; ek (s,Q~V(j(s))^U 



-c 



where s stands for the arc-length parameter of 7, and ( for a system geodesic coordinates normal to 7. 

Corollary 1 1 . 3| gives also some criterion for the applicability of Theorem |1.1| in fact, starting from a non- 
degenerate geodesic in the weighted metric, via the implicit function theorem for A sufficiently small one 



obtains a curve for which (15 1 and the invertibility of (36 1 hold. In particular, when V is constant, one 



can start with non-degenerate close geodesies on M in the ordinary sense. 
The full proof of Theorem |1.1| will be given in the second part of the paper, [40 . In this first part we derive 



some formal expansions of equation (NLS £ I for ip of the form ( 10 1, and check that the assumptions of the 



theorem (the stationarity and the non-degeneracy conditions) guarantee to find approximate solutions 
up to any powers of e. At a formal level, we consider expansions with coefficients depending smoothly on 
the variable s, and the only obstructions to an iterative solvability of ( NLS £ ) are given by the presence 
of a kernel in the linearization of (|8|, which is generated by the functions (d Xj U), j = 1, ... ,n — 1 
and by ill: this kernel arises naturally from the invariance of (J8j) by translation in K™ and by complex 
rotation. However, we can guarantee approximate solvability up to any order provided 7 is stationary 
and non-degenerate: precisely, we prove here the following weaker version of the above theorem. 

Theorem 1.4 Suppose the assumptions of Theorem \l . l\ hold. Then for any m € N and for e small, there 
exists a function ipe,m with the profile (10 1 and which solves ^ up to order o(e m ). 



6 



As discussed before, also some fast-oscillating functions (along 7) contribute to generate some resonance, 
but we do not discuss this aspect here. Postponing the description of the rigorous proof to the introduction 
of [40], here we limit ourselves to mention the main new difficulty compared to the results in [24], [39] 
[412] and [33] . In our case the solutions are complex- valued, and this causes an extra degeneracy in 
the equation, due to its invariance under multiplication by a phase factor. As a consequence, we have 
a further (infinite-dimensional) approximate kernel, corresponding roughly to factor of ip £ in the form 
e — »/i(s) ; f or j 1 arbitrary. The correction in the phase can also be determined by a formal expansion in 
powers of e and, as for /', we still obtain nonlocal terms. Also, when expanding formally the solutions 
in powers of e, the highly oscillatory behavior of solutions generates an increasing number of derivatives 
in s: anyway in the first part, where formal expansions are carried out, this aspect is not very relevant. 

The plan of the paper is the following. In Scction[2]we study the functional in ( 14 1 constrained to the class 
of curves T, and we determine the Euler-Lagrange equation together with the non-degeneracy condition. 
In Section [3] approximate solutions to (NLS e I up to order e are found, and the error terms of order s 2 are 



displayed. The functions we obtain are allowed to depend on a section $ of the normal bundle to 7 and 
on a scalar function fi . These correspond to some tilting of the approximate solution perpendicularly to 
the limit set and to a variation of the phase, in order to have more flexibility. In Section [4] we consider the 
terms of order e 2 and we choose /1 and $ so that even the terms of order e 2 in the expansion vanish: we 
then arrive to the proof of Theorem |1.4[ Finally in Section [5] we collect some technical material, namely 



some integral identities and the verification of ( 36 1 



The results in this paper and in (40] . together with the main ideas of the proofs, are briefly summarized 
in the note [41 j . 

Notation and conventions 

Dealing with coordinates, capital letters like A, B, . . . will vary between 1 and n while indices like j, /, . . . 
will run between 2 and n. The symbol i will always stand for the imaginary unit. 

For summations, we use the standard convention of summing terms where repeated indices appear. 

We will choose coordinates (#1, • • • ,x n ) near a curve 7 and we will parameterize 7 by arc- length letting 
x\ = s. Its dilation 7 e :— ^7 will be parameterized by s — ^s. The length of 7 is denoted by L. 

For simplicity, a constant C is allowed to vary from one formula to another, also within the same line. 

For a real positive variable r and an integer m, 0(r m ) (resp. o(r m )) will denote a complex-valued quantity 
for which - remains bounded (resp. ^ tends to zero) when r tends to zero. We might also 

write o e (l) for a quantity which tends to zero as e tends to zero. 



2 Study of the reduced functional 



In this section we consider the functional in the right-hand side of ( 14 1 defined on the set T, representing 
the approximate energy E e of a function concentrated near 7 with the profile (10 1. We first introduce a 
convenient set of coordinates near an arbitrary (smooth) closed curve in M. Then, using these coordinates 
we write the Euler equation and the second variation formula at a stationary point. 

2.1 Geometric preliminaries 

In this Subsection we discuss some preliminary geometric facts, referring for example to |25j . |53j . Given 
an arbitrary simple closed curve 7 in M, we choose coordinates x\ . . . , x n near 7, called Fermi coordinates 
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in the following way. Wc let x\ parameterize the curve 7 by arc-length. At some point q in the image of 
7 we consider an orthonormal (n — l)-tuple (Y2, ■ ■ . , Y n ) which form a basis for Nqj, the normal bundle 
of 7 at q. We extend the Yj's as vector fields along 7 via parallel transport along the curve with respect 



N , namely by the condition 1 
Next we parameterize a point near 7 using the following coordinates (s, y) € R x 



to the normal connection V , namely by the condition Y; = for I — 2, . . . , n 



(s, 2/2, ■ ■ • , 2/n) >-> exp 7(?) (y 2 F 2 H 1- 

where exp 9 is the exponential map in M through the point q. In this way, fixing s, each curve i 1— > iy, 
for y G \ {0} and t close to zero, is mapped into a geodesic in M passing through 7(5). 

Let us now define the vector fields E\ = and J5j = ^ for i = 2, . . . , n. We notice that on 7 each Ei 
coincides with Yi, while E\ on 7 is nothing but 7. By our choice of coordinates it follows that V eE = on 
7 for any vector field E which is a linear combination (with coefficients depending only on s) of the Ej's, 
j = 2,.. In particular, for any l,j — 2,...,n, and for any aelwe have Vei+uE- {Ei + aEj) = 
on 7, which implies V E t Ej + V e,Ei — for every l,j = 2, . . . ,n. Using the fact that Eas are coordinate 
vectors for A — 1, . . . , n and in particular V e a Eb = V e b Ea for all A, B — 1, . . . , n, we obtain that 
V E t Ej = for every l,j = 2, . . . , n. This immediately yields 

$n$0 = E m (Ei 1 Ej) = (V E m Ei, Ej) + (E u VE m Ej) = on 7, Z,j,m = 2 } ...,n. 

Moreover, still since the -EU's are coordinate vectors for A — 1, ... ,n, we obtain 

d m gij = E m {E 1 ,Ej) = {Vs^uEj) + (Ei,V Em Ej) 

= (V El E m ,E J ) + (E 1 ,V Em E 3 ) = on 7, m,j = 2,...,n. 

Here we used the fact that E m = on 7, namely that Ve 1 E„ 1 has zero normal components. 

If H = H m E m is the curvature vector of 7 (which is normal to the curve), then one has (V^-E^, E\) — 
—H rn on 7, so we easily deduce that 

(16) d m gu = E m (E u E 1 )=2(V El E m ,E 1 ) = -2H m on 7. 
One can also prove that the components Ri m ij of the curvature tensor are given by 

(17) R lmlj = -\d 2 m9ll + H m W. 



Indeed, we have 



-Rimij = (R(Ei, Ej)E\, E m ) — (Vej VejEi, E m ) — (V^ Ve x Ei, E„ 

= (Vbj V Ej Ei , E m ) — Ej (Vei Ei , E m ) — (Vsi E\ , Ve, E m ) 

= (Vfij V ' Ej Ei , E m ) — Ej (Vbj -Ei , £ m ) 

= (Vbi ^ EjEi, E m ) — EjE 1 (Ei, E rn ) + Ej(Ei,\7 El E m ) 

= {Ve^e^E^+E^E^Ve^) 

= E 1 (X7 Ej E 1 , E m )-(X7 Ej E 1 ,X7 El E m ) + ^E J E m {E x , E l ) 



where here we have used the above properties and the fact that 

VfyEh. = V E ,Ej = -d^xxEi. 



Using (16 1 and (17), the above discussion can be summarized in the following result. 



8 



Lemma 2.1 In the coordinates (s,y), for y close to zero the metric coefficients satisfy 

n 1 n 

1 - 2 ]T H m y m + 2 E ( H ' nW - R irnijU)y m yi + 0(\y\ 3 ); 

m=2 m,l=2 

1 n 

, E d ™i9ij\iV™yi + 0(\y\ 3 ); 



m{y) 
9ij(y) 
9kj(y) 



mj.=1 

Skj 



1 " 

J E ^ 

i9kjUy m yi + 0(\y\ ) 



,1=2 



The second derivatives d^gij and df nl gkj could be expressed in terms of the curvature tensor and the 
curvature of 7 reasoning as for ( 17 1. However for our purposes it is not necessary to have such a formula, 
so we leave the expansion of these coefficients in a generic form. 

2.2 First and second variations of the length functional 

We recall next the formulas for the variations of the length of a curve with respect to normal displacements. 
We start with a regular closed curve 7 in M of length L, which we parameterize by arc- length, using a 
parameter s £ [0,L]. Then we consider a two-parameter family of closed curves 7* li t 2 : [0,L] — * M, for 
ti,i2 hi a neighborhood of in K, such that 70,0 = 7- The length L{t\,t2) of 7t ll t 2 is given by 



L(h,t 2 ) = 



dl = 



7*, ,i 2 



(jt 1 ,t 2 ,'jt 1 ,t 2 ) 2 ds, 



where dl is the arc-length parameter and 7t 1> t 2 stands for d7 ^L'* 2 . We also define the vector fields V, W 

along 7t 1) t 2 as V = a T^' 2 and W = ^T^' 2 ■ In the above coordinates, the vector fields V and W along 7 
can be written as 



v = Ev j '(5)^; 

i=2 



W = E W™(5)^r, 



m=2 



(18) 



Differentiating L(t\ : t2) with respect to *i we find 

dL{t 1: t 2 ) _ /"^ (Vy7*i,t 2 ,7*i,* 2 > 



n (7*i,* : n7*i,* 2 ) 2 



Using (16 1, at (ii,^) = (0,0) we have 

(Vv7t 1 ,* a ,7tx,t 3 >=-V m fT n , 
therefore we can write the variation of the length at 7 in the following way 
dL(t u t 2 ) 



(19) 



dtx 



l(*l,* 2 ) = (0,0) 



V m H m ds 



(v,n)ds. 



Using (|T8| we can evaluate the second variation of the length as 

(Vw7t! ,t 2 ) v v7ti ,t 2 > + (7*i ,t 2 > v w V v7ti ,t 2 ) (7*i ,* 2 > v v7*i ,t 2 > (7*1 ,t 2 > v W7*i ,t 2 ) 



d 2 L(h,t 2 ) 
dttdh 

so at (ti 7 t 2 ) = (0,0) we find 
a 2 L(ii,t 2 ) 



(7*i,* 2 ,7*i,t 2 ) = 



(7*i,* 2 »7*i,* 2 ) i 



<is, 



l(*i,t 2 ) = (0,0) 



[( Vw7, V V 7> + (7, V w Vv7> ~ (7 j v v7) (7, Vw7>] ds. 
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Using the definition of the Riemann tensor and the fact that V and W are coordinate vector fields (so 
that [V, W] = 0) the last formula yields 



d 2 L(t u t 2 ) 
dhdt2 



l(ti,t 2 ) = (0,0) 



[(V^W, V) + (7, V W V 7 V> - < 7 , V) (7, V-yW)] ds 



[(V^W, V) + (R(W, 7) V, 7) ~ <7, V^V) (7, V~W>] 



/ (V-yy, V w V)ds. 

JO 



Here, we have used the fact that 5(7, VyyV^V) = (i?(>V, 7)V, 7) + 7<V W V,7) - (V^7, V W V) and 



Jo' 7(^w^> 7)ds = 0. Since V^Ej = on 7 for Z, j = 2, . . . , n, we have 



(Vt,7,V w V)c?s = 0. 



Moreover, recalling (16) we obtain 



V V = ^ V% + X V'Vr./-/, = ^ V'/-:, - X //'V/v 

i=2 j=2 j=2 j=2 



This implies, at 7 



(V^V, V^V) + (R(W,j)V,j) - (7, V^)( 7 , V^W) 
In this way the second variation of the length at 7 becomes 



i=2 j,l=2 



(20) 



dtxdt 2 



l(ti,t 2 ) = (0,0) 



„L / " n \ 

/ - 51 ftijnV J 'W* da. 

7o V i W=a / 



2.3 Determining the phase factor 



In this section we derive formally the asymptotic profile of the solutions to (NLS E \ which concentrate 
near some curve 7, and we determine some necessary conditions satisfied by the limit curve. For doing 
this, using the coordinates (s, y) introduced in Subsection 2.1 we look for approximate solutions ip(s,y) 

'k{s)y s 



of ( NLS e \ making the following ansatz 



<l>{s,y)=e- l eh{s)U 



3€[0,£], ye 



where the function U is the unique radial solution (see [15], [27], [34], [54]) of the problem 
(21) 



-AU + U = U p inR"- 1 ; 
U(y)^0 asM^+oo; 



U > in R" , 

and where the functions /, h and k are periodic on [0, L] and have to be determined. With some easy 
computations we obtain 



d± = if (s) 
ds e 



h(s)U 



k(s)y 



e e + e c h (s)U 



k(s)y 



+ e- l zh{s)k'{s)V v U 



H s )v \ V. 
e ) e ' 
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d 2 1p 



k{s)y 
e 

k(s)y 



-2i f ^h'(s)U 

£ 



k(s)y\ y 



h(s)(k'(s)) 2 V 2 y U 



e 

k(s)y 

e 



+ 2h'{s)k'{s)WU 
+ h"{s)U 



k(s)y\ y 



y y 

e e 



£ 

k(s)y 

£ 



+ h(s)k"(s)V v U 



k(s)y\ y 

e 



and also 



A y ip(s,y) 



(*(3)) a 



■AM 



k(s)y 



-h(s). 



Since U decays to zero at infinity (exponentially indeed, by the results in |27p. and since the function 
-0 is scaled of order e near the curve 7, in a first approximation we can assume the metric g of M to 
be flat in the coordinates (s, y), see the expansions in Lemma 2.1 We look now at the leading terms in 
(NLS £ I, which are of order 1. Since —A g i[i is multiplied by £ 2 , we have to focus on the terms of order ^ 

in the Laplacian of if;. In the above expressions of ^=r and A y ip, we have that the function U and 
its derivatives are of order 1 when \y\ — 0(e), therefore when the variables y appear as factors in these 
expressions, we consider them to be of order e. For example, V 2 U ^ k ( s J y ^ [i^ w jH be regarded as a 
term of order 1. 

With these criteria, using the above computations and assumptions, imposing the leading terms in 
( NLS S \ to vanish we obtain 



-k 2 (s)h(s)A y U 



k(s)y 



+ h(s) [V(s) + (f'(s)) 2 ]U 



k(s)y 



h(s) p U 



k(s)y 



From (21 1, we have the two relations 
(22) k 2 (s)=h{s) p - 1 ; 



[V(s) + (f'(s)) 2 ]=k(s) 2 = h(sy-\ 



We next obtain an equation for /, which is derived looking at the next-order expansion of (NLS e ). The 
next coefficient arises from the terms of order j in —A g ip, which are given by 



f"(S) 



h{s)U 



Multiplying this expression by U 
as well gives 

= f"(s)h(s) 



£ 

k{s)y 



-h'(s)U 



-h(s)k'(s)V y U 



k(s)y\ y 

£ 



and integrating in y € 



imposing vanishing of this integral 



U 2 



k(s)y 



dy + 2h'(s)f'(s) / U 



k(s)y 



dy 



+ 2f(s)h(s)k'(s) / U 



k(s)y 



k(s)y\ y 



dy. 



Integrating by parts and reasoning as for the usual Pohozaev's identity we obtain that / must satisfy 

f"(s)h(s) + 2f'(s)h'(s) - (n - l)f'(s)h(s) ^ = 0. 



This is solvable in /'(s) and gives, for an arbitrary constant A 



(23) 



(n-l)(p-l) 



f'(s) = Ak(s) n - L h(s)- Z = Ah(s) — 



where we have used the above equation (22 1 for k. Now we can solve the equation for h(s) depending on 
the potential V(s) and the above constant A. In fact, we get that h(s) should solve 



(24) 



V(s) + A 2 h{s) 2<J := V{s) + A 2 h(s) 



(n-l)(p-l)-4 _ / 1 (s)? ) - 1 
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where we have set 



fl*\ (n-l)(p-l) 
(25) <T= 2. 

Remark 2.2 FFe notice that, assuming A to be small enough (depending on V and p), the above equation 
is always solvable in h(s). More precisely, when p < ^| (and hence when 2er < p— 1), the solution is 
also unique. For p > i/iere might be a second solution. In this case, we just consider the smallest 
one, which stays uniformly bounded (both from above and below) when A is small enough, see Figures [7] 
anrf[H below. 




Figure 1: the graphs of V + A 2 h 2tT and hP~ x for p < 2±| and for p = ^| with A < 1 




Remark 2.3 In £/ie above expansions, considering the terms of order e, as already noticed, we considered 
the metric g to be flat near the curve 7, and we tacitely assumed the potential V to depend only on the 
variable s. Indeed, expanding the Laplace- Beltrami operator and the potential V taking the variables y 
into account, we obtain an extra term of order e which does not affect our computations since it turns out 
to be odd in y, so it vanishes once multiplied by U ^Mfl^ an( f integrated over K™ -1 . For more details, 
we refer to Sections where precise estimates are worked out (in a system of coordinates scaled in e). 
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2.4 The Euler equation 

Using the heuristic considerations of the previous subsection, we now compute the energy of an ap- 
proximate solution ip concentrated near a closed curve 7, and then find the 7's for which this energy is 
stationary. We let ip 7tA denote the function constructed in Subsection 2.3 In order for the function 7i ^t 
to be globally well defined, we need to impose one more restriction, namely that is periodic in s 

with period L. This is equivalent to require that f'(s)ds is an integer multiple of 2ire, since we have 



the phase factor e 1 « 
integer multiple of 27re. 



Multiplying §NLS e 
solution we find 



in the expression of "07,-4- From (23 1, then we find that also J Q h(s) a ds is an 
by 07,-4 and integrating by parts, from the fact that ^7, -4 is an approximate 

1 



M 



( £ 2 |V^ 7 ,^| 2 



1 



p + 1 



A I 



\^, A \ p+1 dV g . 



dV g - 



P 



1 



M 



IV>7,. 



IP+1 



Since V> 7 ,-4 is highly concentrated near 7, using the coordinates (s, y) introduced in 



2.1 



we have that 



M 



\^, A \ P+1 dVg 



ds 



h(S) 



p+1 



U 



k(s)y 



P+i 



dy. 



Using a change of variables, the last two formulas and (22 1 we find that 

E e (^, A ) -Cos 71 - 1 



(26) 
where 

and where we have set 
(27) 



C = 



= P 



1 



h(s) ds, 



\U(y)\ p+1 dy 



l--(p-l)(n-l)=p-a-l. 



Consider now a one-parameter family of closed curves 7 t : [0, L] — > M , where t belongs to a neighbor- 
hood of in K and where 70 = 7. We compute next the approximate value of the derivative in t of the 
corresponding energ y de fined by ( 26 ) 
As in Subsection 



2.2 



we let Vt denotes the vector field Vt(s) = ^Sr (s) and we assume that V := Vo is 



normal to 7. For any t near zero, we let fct(s), ht(s), ft(s) be defined by (22 1 replacing 7 by j t and V(s) 
by Vt(s) :— V(7t(s)). Since we require periodicity of each curve 74 in the variable s, we also allow the 
constant A given in (23 1 to depend on t. Denoting this by At, by the above considerations we choose At 
so that the following condition holds for every value of t 



(28) 



Atht(s) a ds 



ft(s)ds — const. 



Below, we let A' t = ^At and we will consider h t (s) as a function of At while Vt(s) as implicitly defined 



in (24 1. From (19 1 and 

i-L 



dVtjs 
dt 



Acrh- 1 ^(V N V,V)ds-A 



| t=0 = (V Ar U(s), V(s)), differentiating (28 1 with respect to t at t = we get 

L 



h a (V,H)ds + AA'<T / h 



dX 



ds + A' / h a ds = 1 



where we have set A' = A' and where V W V stands for the component of VU normal to 7. From this 
formula we obtain the following expression of A' 



(29) 



A' 



-A 



J a L (ah^§^(V N V,V)-h"(V,H))ds 



Jo (Aah° 



BA 



h a ) ds 
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Similarly, computing the derivative of the (approximate) energy with respect to t we find 



dE e {u ltAt ) 
dt 



eh' 



t-idh 
dV 



W N V,V) -h a (V,B.) +9A'h 



dA 



ds. 



Using ( 29 1 we deduce that the variation is given by 

rL dh 



dE £ (u Jt , At ) , 
dt 



o dV 

L 



V iV U,V) 



0h» 



J L {Aah'-^ + h-)da 



ds 



(V,H) 



Differentiating ( 24 1 with respect to A and V we get 

2Ah 2,T 



(30) 

so it follows that 
(31) 



Ah"J L 9h e -^ds 
f L (A*h"-^ + h")ds 

_ = = 2Ah 2 ° — 

dA (p - l)hP~ 2 - 2aA 2 h 2 °- 1 <9V 



ds. 



Oh 



Aah^^ + h 
dA 



(p- l)hP~ l 



(p-l)h - 2oA 2 h°' 



Similarly, since 9 = p — a — 1, see (25 1 and (27 1, we deduce that 

2AhP+ a - 2 



Therefore we find also 
(32) 



_! dh _ 

dA ~ (p-l)hP- 2 -2aA 2 h 2 °- 1 
2A9 



Oh6-l dh 
011 OA 

Acrh—^ + h 



p-1 



Hence from the last formulas the variation of the energy becomes 



(33) dE ^^\ t= , 



dt 



dV 



Oh" - 2 -^h^ 



p-l 



ds 



/ L (V,H> 
Jo 



2 A 2 9 . 
p~^l' 



ds. 



Also, from the second equality in (30), dividing by h a , multiplying by and using the identity p—a— 2 
9 — 1 we obtain 



p-l 



2aA 2 h^ dh 



dh 



_i dh 



dVp-l 



dV dV 



Using ( 33 1 and the last formula, we get the following simplified expression 

rL a 



dt 



\t=o 



o P- 1 



V^V.V) - (V,H) 



^-Ah"- 1 - 2A 2 h 2 ° 



ds. 



Therefore the stationarity condition for the energy (under the constraint (28 1) becomes (V V, V) 
(V, H) (Jt^hP^ 1 — 2A 2 h 2a ) for every normal section V, namely 



(34) 



V N V = H 



p-l 



h p - 



2A z h 



2,2a 



We will see that this formula will be crucial later on to find approximate solutions. 

[p- l)hP~ l 



Remark 2.1 By (31 1, we have that 



d 

dA = (p - l)h e - 2oA 2 h° ' 

If A is sufficiently small (depending on V andp), then we have A (Ah a ) > 0. This will be used in the 
second part \40^ when, for a fixed e, we will adjust the value of the constant A for obtaining periodicity 
of the function f . 
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2.5 Second variation and non-degeneracy condition 

We evaluate next the second variation of the Euler functional. As in Subsection 12.21 we consider a 
two-parameter family of closed curves 7i ll t 2 , where £i,£2 are two real numbers belonging to a small 
neighborhood of zero in R, and where 70,0 = 7- As before, we require the constraint (28 I along the whole 
two-dimensional family of curves, and we assume the functions f,h,k and the constant A to depend on 
£1 and £2, an d we will use the notation At lt t 2 , etc.. Keeping this in mind, we define the two vector fields 
^*i>*2 = ^dti 2 ' Wfr,,t 2 = a ^ t2 1 and we can assume that V := Vo,o, W := Wo,o are normal to the initial 



9*2 



curve 7. With some computations, which are worked out in Section 5.2 one finds that, at (£i,£2) = (0,0) 
d 2 E £ (u^ ut2 , Atu J = r L L 2A 2 6 h 



<9£i<9£ 2 



ds 



(35) 



ad 
P- 1 Jo 



- l ^{V N V,V){V N V,W)ds 



29 
P~ 1 Jo 



Here, A\ stands for ^%f 3 -|(t 1 ,t 2 )=(o,o), V J = and = where the V J 's, W J 's are the compo- 

nents of V and W with respect to the basis (Ej)j introduced in Subsection 



ds. 



2.1 



Integrating by parts and using ( 29 1 , from the last formula one derives that the non-degeneracy con- 
dition is equivalent to the invertibility of the linear operator Z ■ x(-^l) ~~ * xC^7) (from the family of 
smooth sections of the normal bundle to 7 into itself) whose components are defined by 



(W 



p- 1 



v 



- e h - 1 - 



2A 2 a 
/ 

p-l 



+ -h-° {((V") V)(V, E m ) - H m (V N V, V) - (H, V)(V N V, E m )} 

p — 1 

- -^rh-W^-(V N V, V)(V N V,E m ) —AA\ (vh°-^(V N V,E m ) - h°H>») , 
p — 1 ov p — 1 \ / 

dh J^ ■ Recalling formula (17), using (34) and some other elementary computations we obtain 




V m - 9 h - 



2A 2 a.. 
p-l' 



h"- 1 tiv rn + — -r'((v w )V)(v,£,„) 



+iJ m (H,V) 



26» 
~P~^l 

2A 2a 



AA\ ( ah 



p-l 

.idh 
dV 



(V N V,E m )-h°H r ' 



(p-l) (3+ f)h e -8aA 2 h a 
(p- l)h e -2A 2 ahF 



For future convenience, we expand the last product explicitly, finding 

'-(p - 1) (3 + f ) /V 29 - l^M^h 2 " + 2A 2 (5a + 39)h e +° 



iJ m (H,V) 
We also notice that 

ah a - 1 (V N V,E n 



(p-l)h e - 2A 2 ah° 



dh (p-\)(9-a)h^ 
8V 9[(p-l)h e -2aA 2 h°] 



H' 
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In conclusion the non-dcgcncracy condition is equivalent to the invertibility of the operator 3 : x(-^7) 
x{N"f) given in components by 



(36) 



2A 2 e j 
p - 1 



2 V P-l 



p-l 

E(4n5n)V J I -2.4.4', 



p-l 



h-°((v»yv)[v,E m ] 



H m (H,V) 



-(p-l)(3 + f)/» 



29 _ 16ctMI/,2ct 
p-l 



[{p-l)h e -2a A 2 h°] 
2A 2 (5a + 38)h 8+c 



(p - l)h e - 2A 2 oh° 
We summarize the results of this section in the following Proposition. 



Proposition 2.4 Consider the functional on curves J h (s) ds restricted to the setT in (13 1. Then the 



of the operator 3 in ( 36 1 



stationarity condition is (34 1 and the non-degeneracy of a critical point is equivalent to the invertibility 



3 Approximate solutions 



Using a change of variables, equation ( NLS £ \ is equivalent to the following 



(37) 



Ag E ip + V{ex)ip 



in M F 



where M e denotes the manifold M endowed with the scaled metric g e = j^g. With an abuse of notation 
we will often denote it through the scaling M e = |M, and if x € M e we write ex to indicate the 
corresponding point on M. 



In this section we find a family of approximate solutions to the scaled equation ( 37 1 . We consider a 



simple closed curve 7 which is stationary within the class T, namely satisfying (34 1. First, we introduce 
some convenient coordinates near the scaled curve 7 6 = -7, expanding the Laplace-Beltrami operator with 
respect to the scaled metric in powers of e. Then, using these expansions, we construct the approximate 
solutions solving formally (37) up to order e. Since in the second part |40j we will need to work out 



rigorous estimates, in order not to repeat later the expansions we will treat some terms carefully and not 
only at a formal level. 



3.1 Choice of coordinates in M £ and expansion of the metric coefficients 

Using the coordinates (s, y) of Section [2] defined near 7, for some smooth normal section ^(s) in Nj, we 
define the following new coordinates (s,z) (here and below we use the notation s = es) near -7 



(38) 



y-<f>(ss); 



z e 



In this choice we are motivated by the fact that in general we allow the approximate solutions to be tilted 
normally to j £ , where the tilting $ depends (slowly) on the variable s: this allows some extra flexibility 
in the construction, as in [23], [33] and [32] . As we will see, the choice of $ is irrelevant for solving 



(37 1 up to order e; on the other hand, the non-degeneracy assumption will be necessary to guarantee 



solvability of the equation up to higher orders. 

We denote by cjab the metric coefficients in the new coordinates (s,z). Since y = z + <fr(es), it follows 



9cd 



( d VA\ ( dy B \ 
fej fcj 



Explicitly, we then find 



gu = 9uU+9 + 2ey^$jgijl z+ 3. + e 2 5>;.( £s )<i>' m Ms 

j 
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9lj — 9lj\z+<S> + ®'m( £S )9jrn\z+<S>] 9jni — 9jm\z+<S>- 

m 

At this point, it is convenient to introduce some notation. For a positive integer q, we denote by R q (z), 
R q (z, <&) and R q (z, <&, <&') error terms which satisfies respectively the following bounds, for some positive 
constants C and d 

\R q {z)\<Ce«{l + \z\% 

\R q (z^)\<CrJ(\ + \z\ d ); 
\R q {z^) - R g (z,$)\ < Cs"(l + \z\ d )[\<i> - $|], 

and 

\R q (z,$,&)\<Cei(l+\z\ d ); 

\R q (z, - R q (z, l>, &)\ < Ce«(l + |$ - $| + |$' - <£'|]. 

We also introduce error terms involving also second derivatives of i? g (.z, <&") which satisfy 

|i? g (z, $, $")| < Ce«(l + \z\ d ) + Ce q+1 (l + \z\ d )\<S>"\; 

\R q (z, - R q (z, $")| < Ce 9 (l + |z| d )[|$ - $| + |$' - *'|] (l + E (|$"| + |*"|)) 

+ Ce' 3+1 (l + |z| £i )|$"-l'"|. 



2 

m,Z=2 



Using the expansion of the metric coefficients gAB in Lemma 2 . 1 1 and this notation, we then obtain 

n ^ n 

ffii = 1 -2e ^] i? m (z m + $,„) + -e 2 d 2 ml g ll \ 1 {z m + <^ m ){zi + ^i) 

m=2 

(39) + e 2 \<S>'\ 2 + R z {z, $,$'); 



1 " 



2 

m,Z=2 
1 - 

ffiy = ^fej + -e 2 51 d5fcil 7 («m + $ m)(«J + *i)+ -Ra(*. 

m,Z=2 

Next we compute the inverse metric coefficients. Recall that, given a formal expansion of a matrix as 
M = 1 + eA + e 2 B, we have 

Af~ 1 = 1 - eA + e 2 A 2 - e 2 B. 
In our specific case the matrix A is the following 

(40) A = 

and the elements of its square are given by 

(A 2 ) n = 4 ( £ ir™(z m + <f> m )) + (^) 2 ; 



-2 £ ff m (z m + $ m ) 

m=2 





\m=2 



(A 2 ) l3 = -2 £ iT»(z m + 4> m ) (*;•) ; and (A 2 ) y = (${) 
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Therefore, using the above formula, the inverse coefficients are 

n - n 

= 1 + 2e H m (z m + * m ) - -e 2 ^ ^ i5 n| 7 (2 m + $ m )(zi + $i) 



3 11 



m,i=2 



4e 2 ( £ iT m (* ro + $ w )j +i?3(z, $,$'). 



Vm=2 



We also get 



1 ™ / ™ \ 

m,/=2 \m=2 / 

+ Ra(z, $,$')• 

Moreover 

J=2 

Similarly, with some simple calculations one also finds 

n n 

ft(s n ) = 2e 2 ^(ff™)'(z m + a> m ) + 2 £ 2 ^ff m ^ + i? 3 (z,a>,ci>',a>"). 

m=2 m — 2 

Differentiating now § J with respect to the first variable we obtain 

Vm=2 / 

Analogously, we get 

1 " 



,.1=2 



ft(ff fcj ) - ~e 2 d 2 kl 9kj\M + *i) + ^ *. *')■ 



1=2 



Finally, using the formal expansion gcD = 3cd + £Ac<d + £ 2 Bcd + °(£ 2 ), analyzing carefully the error 
terms we obtain 

Vtetj = 1 + \etr(A) + e 2 (^-{tr{A)) 2 - \tr(A 2 )^ + \e 2 tr{B) + 0(e 3 ). 

From the above expressions in p9j ), ( |40| ) we deduce that 
v/det! = l-e^iT ro (z m + $ m ) 



m,Z \m=2 / 
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d my ^M^ = -eH m + e 2 
+ R 3 (z, $,$'), 

moreover 

m m 

The Laplacian of a smooth function u in coordinates (s, z) has the following expression 

-A~ g u=-Y,~g AB d 2 AB u-J2Mg AB )dAU-^=J2~9 AB (d B V^) d A u. 
a,b a.b y act g ^ g v / 

We are going to expand next each of these terms. First, we consider the determinant of g. Recall that 
for a matrix of the form 1 + eA + e 2 B the square root of the determinant admits the formal expansion 

(41) ^drtg = 1 + \trA + e 2 Q(ir A) 2 - ^tr(A 2 ) + ^trb^j + o(e 2 ). 

Lemma 3.1 Let u be a smooth function. Then in the above coordinates (s,z) we have that 
A- g u = d 2 ss u + A z u-eJ2 H'd.u - 2e ]T ^d%u + 2e{H, z + <S>)d 2 s u 

3 3 

- e 2 (H, z + *) Wd jU - l -e 2 d 2 ml9ll {z m + $ m )(z ; + + 4e 2 (H, z + <£) 2 d 2 s u 

m,j 

- e 2 d 2 ml9lj (z m + <J> m )(z ; + ^)d%u - 4 £ 2 (H, z + $) ^ $^.« + e 2 X] 

- ^ 2 Z! + + *')9 t > + e 2 (H', z + $)9 s u - e 2 Z d %9^ z i + 

m,l l,j 

+ i?g(z, + ^3(«, &)d%u + R 3 (z, $, + i2s(«, *") (d s u + 0,-u) . 

Moreover, given two smooth normal sections $ and $ and defining the corresponding coordinates 

(s,y — $(ss)) and (s,y — &(ss)) 
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and set w$(s, y) :— u(s, y — $(es)), u^(s, y) := u(s, y — $(es)). We then have 
A- g u* - A- g u$ = -2e - & 3 )d%u + 2e(H, $ - $)d 2 ss u + e 2 ~ 



e 2 X! 



■in j. 



-e 2 V a 2 



2z m ($ ; - $j) + $i($ m - $ m ) + $j($ m - l> m ) 
2z m ($ ; - $,) + $ ; ($ m - $ m ) + <^($ ro - $ m ) 



3 S > 



in.! 



2s 2 J2 HJ 



x9u 



z,(${-$0 + *i(*{-*{) + *i(*i-*0 



4e 2 ^ 



in. I 



j t,j,i j 

e 2 (H', $ - $)d s u -e 2 Y, dfjgij^i - 

1,3 



~ 2e 3 H ™ 

mj 

+ 0(l + \z\ d ) 
+ 0(l + \z\ d ) 



5 4 (|$ - *| + |$' - $'|)|a s 2 s u| + £ 3 (|<i> - $| + |$' - &\)Q%jU\ + \d 2 u)\ 



e»(|$ - $| + |$' _ $'|) + e 4 (|$"||$ - $| + |$" - $"|) 



(|cU| + \dju\). 



Proof. The proof is based on the Taylor expansion of the metric coefficients given above. We recall 
that the Laplace-Beltrami operator is given by 

A- g = J2 M Tdet^ (g e ) AB d B ) , 

where indices A and B run between 1 and n. We can also write 

A s = E (~9 AB 9 ab + (9a ~9 AB ) d B + -j^~9 AB (d B y/tet§) B 



A.B 



Using the expansion of the metric coefficients determined above and (411, one can easily prove that 



J2~9 AB d 2 AB u = A z u + d 2 ss u~2eY / ^d 2 3 u + 2e(n,z + <P)d 2 s u + e 2 Y / ^d 2 u 

AB j l,j 

+ 4e 2 (n,z + ^) 2 d 2 ss u--e 2 ^2d 2 nl g kj {z m + ^ m ){z l + ^ l )dl j u 

m,l 

~ \ £2 E d ™i9u( z ™ + *m)(*i + ®i)3%u - 4e 2 (H, z + $)^-9 s > 

m,l 

~ E Oll(*m + $m){zi + + M*, *, (B^U + d%U + df jU ) 
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d A ~9 AB d B u = -e 2 ^d 3 u ~ e 2 £ 8fag kj (zt + $i)d jU ~ 2e 3 (H, z + $) £ 

+ 2e 2 (H', z + $)d s u ~e 2 J2 %j9ij{*\ + + ^(z, + S/tt) 

^-^=g AB [d B \fd^tpjd A u = -eY,H j d J u-e 2 (B.,z + ^)Y^H j d 3 u-e 2 (H',z + ^>)d s v 



A.B 



The result then follows by collecting these three terms. 



3.2 Expansion at first order in e 



In this subsection we solve formally equation (37 1 up to order e, discarding the terms which turn out to 
be of order e 2 and higher. 



For the approximate solution as in (|10|), we make a more precise ansatz of the following form 

s G [0,27r],i/G K™-\ 



(42) ^i, e (s,z) = e 1 e {/i(es){7 (fc(es)z) + e [iu r + no*]} , 
where foiss) = f(es) + efi(es). By direct computation, the first and second derivatives of ipi l£ satisfy 



d s ip 



l.e 



-ifo(es)h(es)U(k(es)z) + eti(es)U(k(es)z) + eh(es)k'(es)VU(k(es)z) ■ z 



r* « He/V + e/V] + 0(e 2 ); 



,/o(= 



[/i(es)A;(es)9i?7 (fc(es)z) + ediW r + iediWi 



-{f' fhU(kz)e~ lM ^ 1 - iee- %M ^ [f'hU(kz) + 2f'tiU(kz) + 2fhk'\7U(kz) ■ z] 



Q 2 * 



ef rz e' lM ^[w r + iw l ] + 0{e 2 ); 



,- /o(") 



[/i(es)fc 2 (es)3; 2 J /7 (&(es).z) + edfjW r + iedfjWi] ; 



.- /o(") 



i/o(es)/i(es) + eh'(es) 



k{es)djU{kz) 



:/0( = 



' h(es)k' (es) 



kJ2dfjU(kz) Zl + djU(kz) 



; J0(») 



;/0i££) 



— ief'(es)djW r (ss, z)e ' + ef (es)d J w i (es 1 z)e ' +0(e 2 ). 
Similarly, the potential satisfies 

F(ea;) = F(es) + f(V N y, z + $) + 0(e 2 ). 



Expanding (37 1 in powers of e, we obtain 

■■/<]("■ 



A ff Vi, £ + V(ea;)Vi ie - |V>i, E | p ~Vi, E ) = eK r + teU, + 0(e 2 ), 
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with 

(43) K r = C r w r + 2f'f[hU + 2/' 2 M7(/cz)(H, z + $) + hk{H, VU(kz)) + (V N V 1 z + ^)hU(kz); 



(44) 



Ki = C t w l + [fhU(kz) + 2f'h'U{kz) + 2f'hk'VU(kz) ■ z] -2^y& j fhkd j U{ki 



and where we have defined the two operators C r and Ci as 

C r w = -A z w + (V + f 2 )w - ph p - 1 U(kz) p - 1 w; 

CiW = ~A z w + {V + f 2 )w - h p - 1 U(kz) p ~ 1 w. 

It is well-known, see for example [51], that the kernel of C r is generated by the functions d2U(k-), . . . , d n U (A;-), 
while that of Ci is one-dimensional and generated by U(k-). 

We choose the functions w r and Wi in such a way that lZ r and IZi vanish. Since C r is Fredholm, 
the solvability condition for w r is that the right-hand side of this equation is orthogonal in L 2 (]R™ -1 ) 
to d2U(k-), . . . , d n U(k-). Therefore, to get solvability, we should multiply the right-hand side by each of 
these functions and get 0. The same holds true for Wi, but replacing the functions d Zj U(k-) by U(k-). 
We discuss the solvability in Wi first. Writing this equation as CiWi — f, we can multiply it by U(k-) and 
use the self-adjointness of Ci to get 



= 



w l C i U{k-) = 



U{k-)CiW. t = 



fU(h). 



Following the computations of Subsection |2.3[ this condition yields 

f"hk- [n - l) + 2f'h'k- {n -^ = (n - l)hk'k- n f, 

which implies 

h.n— 1 



This equation is nothing but (23 1, and hence the solvability is guaranteed. Since Ci clearly preserves the 
parity in z, we can decompose Wi in its even and odd parts as 

Wi = Wi^ e + W^ , 

with Wi te and Wi tQ solving respectively the equations 

CiW he = - [f'hU(kz) + 2f'h'U(kz) + 2f'hk'VU(kz) ■ z] ; 



C i w i , = 2j2[^' j f'hkd j U{kz)] , 



where the right-hand sides are respectively the even and odd parts of the datum in (44 1. We notice that, 
since the kernel of Ci consists of even functions, only the even part of the equation plays a role in the 
solvability, since the product with the odd part vanishes by oddness. 



Indeed, (43) and (44 1 can be solved explicitly, and the solutions are given by 

Wi,o = -J2$' j f'hz j U(kz). 



P —±fh<\z\ 2 U{kz); 



In fact, as one can easily check, we have the following relations 

C l {z j U{kz)) = -2kdjU{kz)\ C % {\z\ 2 U{kz)) = -2{n - l)U(kz) - 4fcW(fcz) • z, 



which imply the above claim (here we also used ( 22 1 and some manipulations) . 
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Turning to w r , if we multiply by djU, we integrate by parts and use some scaling, we find that the 
following conditions holds true, for j = 2, . . . , n 

2H j ((f') 2 [ U{z) 2 dz-^--\ \VU(z)\ 2 dz) +{V N V 1 E ] ) [ U{z) 2 dz = 0. 
V Jr- 1 n-1 J M n-i J Jk"- 1 



Using (22 1, we get equivalently 

U{zfdz~ h —- [ \VU{z)\ 2 dz\ + (V N V,Ej) [ U(z) 2 dz = 0, j 



2iP ^ 2 /i 



= 2,...,n. 



From a Pohozaev-type identity (playing with (21 1 and integrating by parts) one finds 



(45) 



\VU(z)\ 2 dz 



(n-l)(p-l) 



(3-n)(p + l) + 2(n-l)y R „- 1 
Using this formula the solvability condition then becomes 



U{z) 2 dz 



(n-l)(p-l) 
26> 



U{z) 2 dz. 



W (p - 1) 



= (V A V,£7,-), i = 2 



which is nothing but the stationary condition ( 34 1 . Therefore, since we are indeed assuming this condition, 
also the solvability for w r is guaranteed. As for uii, we can decompose w r in its even and odd parts as 

where w r , e and w ry0 solve respectively 

£ r ii; r , e = -2ff[hU -2(f) 2 hU(kz)(H,<S>) - (V N V^)hU{kz)- 

C rWro = -2(f) 2 hU(kz)(H,z) - hk^H j d 3 U{kz) - (V N V, z)hU(kz). 

j 

Using the Euler equation, one gets 

C r w r . = -hJ2 H3 [kdjXJ +hP- lV —^z 3 U 



It is also convenient to have the explicit expression of w r . We notice first that 



Hence it follows 



[h(V N V + 2(/') 2 H, *) + 2ff 1 h] U—^—f Uikz) + ^VCZ(fcz) • z 



Using (|34| we finally find 

-p-l 



h*{H,$)+2ff 1 h 



By the above computations we obtain the following result. 



Lemma 3.2 Suppose h(s) and f(s) satisfy (111 and (12) for some A > 0: assume also that the curve 7 
verifies (151. Then there exist two smooth functions w r (s, z),Wi(s, z) for which the terms lZ r and IZi in 
(43 1 -( 44 ) vanish identically. Therefore, the function ipi tS in (42 1 satisfies (37l up to an error 0{e 2 ). 
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3.3 Expansions at second order in e 

Next we compute the terms of order e 2 in the above expression. Adding a correction e 2 [v r + ivi] to the 

function in ( 42 1 we define an approximate solution of the form 

(46) 



-4> 2e (s,z) = e^ 111 ^ 11 {h(es)U (fc(es)z) +e[w r + iwi] + e 2 [v r + ivi]} ; s e [0,2n],y € E" l , 

where /q = /(es) + e/i(es). The first and second derivatives of ^2,e are given by 

-ifo( L es)h(es)U(k(es)z) + eti (es)U{k(es)z) + eh(es)k' \es)VU(k(es)z) ■ z 
-ief' Q w r + efoWi + -ie 2 f' Q v r + e 2 f' Vi + e 2 (d s w r + iwCj + C»(e 3 ); 



;/0l£f) 



djtp2,e = [/i(es)fc(es)9jZ7 (fc(es)z) + edjW r + iedjWi + e 2 djV r + ie 2 djVj\ ; 



ds 2 



= -(fo) 2 hU(kz) -ie fo'hU(kz) + 2f^h'U(kz) + 2fohk"VU(kz)-z -efg[w r + 



e 2 fo[v r + iv t ] + e 2 2f d sWl + ti'U(kz) + 2h'k'VU ■ z + hk"VU ■ z 



+ hk' 2 V 2 U(kz)[z,z} + ftwi 



l£ 



2ftd s w r + ftw r + 0(e 3 ); 



dfjipi^ = [h(es)k 2 (es)df j U (k(es)z) + edf w r + ied 2 j w l + e 2 d 2 j v r + ie 2 d 2 j v i ] ; 



: M££l„2 , 



— ifQ(ss)h(es)k(es) + eh'(es)k(es) + eh(es)k'(es djU(kz) 
h(es)k(es)k'(es) d 2 jU(kz)zt + djU(kz) — ief' djW r 



+ efQ(es)djWi - ie f djV r + e f' (es)djv l + e d sj w r + ie d 



We also have the formal expansion 



, : hPp^U + pehP-^U^Wr + ieh p ~ 1 \U\ p - 1 w l + -p(p- l)e 2 h p - 2 \U\ p - 3 Uw 2 r 
1 



+ -(p - l)e 2 h p ~ 2 \U\ p - 3 Uw 2 + i(p - l)e 2 h p - 2 \U\ p - 3 Uw r Wi 
+ pe 2 h p - 1 \U\ p - 1 v r + ie 2 h p - 1 \U\ p - 1 v i + 0(e 3 ). 
Similarly, expanding V up to order e 2 , we have 

V(ex) = V(es) + e{V N V, z + $) + ^e 2 (W N ) 2 V[z + z + $] + i? 3 (z, $). 



Using the expansions of Subsection |2.3| we obtain 



,;M£f) 



2.5 



= £ 2 (7l r + i7l;) 

= e 2 (R r ^ e + R r ^ ) + e 2 i(Ri^ e + Ri^ ) 

+ e 2 (R rye J 1 + Rr,o,fs) + £ 2 *(-^e,/l + ^i,o,/i) 

+ e 2 C r v r + e 2 id i v i + 0(e 3 ), 
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where 



R r ,e = -\{f l ) 2 hU{kz)Y J dL9ii{z m z l + ^ m ^ l )+2{f l )\^,w r ^ + w r ^ z) 
+ 4(f') 2 hU(kz)[(H,z) 2 + (U^) 2 ] 

- [h"U(kz) + 2h'k'VU{kz) ■ z + hk"VU{kz) ■ z + h{k') 2 V 2 U{kz)[ Zl z]] 



(47) + 2fd sWi , e + f"w i , e + 2^ j fd j w i , 



+ hkJ2 d 2 m g m3Zl d 3 U{kz) + J2 H m d m w r . + hk(U, z) J2 H m d m U{kz) 



hk 2 d?jU{kz) 



+ khJ2 



(H,z)H m - W%nxmz x 



d m U(kz) 



\p{p iw- 2 u{kzy- 2 {wl e + wl a ) - i(p - iw- 2 u{kzy~ 2 {wl e + w 2 Q ) 

+ (V N V, W r , Z + W r , e <S>) + 1^2 d ™3 V ( Z ™ Z 3 + ^m^j)hU(kz); 



m>3 



R r ,o = -(fThU(kz)J2dLgiiZrn<S>i+8(f) 2 hU(kz)(H 7 z)(H 7 <S>) 

l.m 

+ 2(/') 2 (H, Wr>eZ + w r , <5>) - 2f'd s w h0 - f"w h0 + 2tikJ2& j d j U{kz) 



2h'kJ2^ 



k^2dfjU{kz)zi + djU(kz) 



+ 2f &jdjWi,e + hk(U, *) H m d m U(kz) 



(48) + 



\ d ^m9tj { z m®l + Zl$m) 



hk 2 d 2 tj U(kz) + hkJ2®PjU{kz) + d fm9m^i hkdjU{kz) 



+ J2 Hmd ™ W r,e + hkJ2 



d m U(kz) 



- p(p - l)h p - 2 U(kz) p - 2 w r . e w r . - (p - l)h p - 2 U(kz) p - 2 Wi, e Wi, 
+ (V N V, w r . e z + w r , $) + fyVzjfyhUikz); 

3,1 

Ri, e = 2[f"hU(kz)+2ftiU(kz) + 2fhk , VU(kz)- z] (H, *) 

+ 2(/') 2 (H, Wi e § + Wl o z) + 2f'd s w r , e + f"w r , e - 2f ^'jdjW^o 



(49) 



2f'hkJ2d j U{kz) 



l.m 



2(H, z)^'j + d lm9ij{z m $l + zi$ m ) 



- f'hU(kz) (j2(d 2 m gu* m - 2(H, *'))^ - f'h 
+ I f'h (^2 du9n^ U{kz) + ^ W8 jWit0 



U{kz) 



(p - l)h p 2 U(kz) p - 2 (w r , e w lte + Wr^.o) + (V iv V, Wi, Q Z + Wi te <5>); 
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(50) 



2 [f'hU(kz) + 2f'h'U{kz) + 2f'hk'VU(kz) ■ z] (H, z) + ^ 

i 

2(/') 2 (H, w ite z + w^) + 2f'd s w r>0 + f"w r , - 2f Vftw^ 



- 2f'hkJ2d J U(kz) 



2(H, + i ^ E^g^^ZmZi + <S>i<£ m ) 



l.m 



- fhU(kz) dl m g lxZr )j - f'h df ]9l3 z}j U(kz) + l -fh d 2 ll9ll z^j U(kz) 

- (p - l)h p - 2 U(kz) p - 2 (w rie w tt o + Wr t0 w ite ) + (V N V, w i>e z + w it0 <f>). 

We used the notation Rj.^^, i?,.^^, Ri, e ,fi & n d Ri,ojt for the terms involving f\, namely 
Rr,e, h = (/0 2 ^ + 2/7i«V, e + 4(H, - 2p(p - l)^- 2 |C/|f- 2 /i 2 /' 2 /f f> 2 

>-l) 2 



(51) 



2f'f[h(V»V, *)U + 4<H, f[hU - 2p 



-h 2 P- 1 ff[(H,$)U p - 2 U 2 



(52) 



Br.o.A = 2/7^ I%0 + 4(H,z)/7^C/-2p(p-l)/7^ f, - 1 C/ p - 2 C/ Wri0 
+ 2f'f[h(V N V, z)C7 + 2H\ff l hkd J U + 4<H, z){f') 2 f f[hU; 



(53) 



^, e)/l = 2h'f[U + 2hf[k'VU ■z + 2f'f[w^ + f['hU + 4f'd s (hf'f' 1 U) 
+ 2f"hff[U-2(p-l)hP- 1 \Ur 2 ff[Uw he] 



(54) 4^ - 2/7{ Wi , - 2(p - l)/i p - 1 |^| p - a /'/I^ - Kffhkf^'^U - 2hkf[$' j d j U, 

where we wrote for brevity U — ^ p _^j hP -i U(kz) + ^WU(kz) ■ z. Again, we collect the results of this 
section in one proposition. 

Proposition 3.3 Suppose are smooth functions on [0,L], let z be the normal coordinates given in 

(381 and Wi,w r be be as in Lemma 3.2 Then, if ip2,s is as in (46l, in the coordinates (s,z) we have 



(55) 



-Ag e -02,e + V(ex)lp 2 ,e ~ \^2,e\ P l ^2,e = e 2 (£ r V r + iCjUj) 
+ S 2 (R r ^ e + R r ^ + iRi^e + iRi,o + Rr,eJ\ + Rr,o,f 1 + iRi,e,f 1 + iRi,oJ t ) + o(£ 2 ), 



where the above error terms are given respectively in (47 1 -(54). 



4 Final expansions and proof of Theorem 1.4 



In this Section we prove existence of approximate solutions to (NLS £ I up to any power of e. To do this 



we first evaluate the projections of the error terms on the kernels of the operators Ci and C r , and adjust 



<3> and fx so that these projections vanish, in order to obtain solvability of (37) up to order e 2 . Then, 
with an iterative procedure, we turn to the general case. 

4.1 Projection onto the kernel of Li 



If one wants to find Uj so that the imaginary terms in (55 1 vanish, by Fredholm's alternative the imaginary 



R's must be orthogonal for every s to the kernel of Ci, which is given by iU{k{es)-). To compute this 
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projection, by parity reasons, we need to multiply i?i je and i?i, e ,/i by U(k(es)-) and to integrate over 
R n_1 . We evaluate the two terms separately. 



Contribution of i?i >e - After some manipulation we obtain 

Ri,e — C\ + C*2 + C3 + C4 + C5, 

where 

C y = 2 [f'hU(kz) + 2f'h'U(kz) + 2f'hk'VU(kz) ■ z] (H, $) 



C 2 = -2f'hkY,djU{kz) \J2dim9iiZm$i - f'hJ2df l9l ^iU(kz); 



l,m 



j.i 



2(H, $') - \ (j2 d 2 u gu®l) J ni U(kz) 



C 3 = -A^fhkdjUikz)^^)^ - 2fh^H j ^' ] U{kz)- 

3 3 

C 4 = -(f) 2 (H,w l . e ^ + w li0 z)+Y,H l d l w l . 

1 

+ (V N V,w i>e $+w ii0 z) - (p- l)h p ~ 2 U(kz) p ~ 2 (w r>e w iie + w r , w it0 ); 
C 5 = 2f'd s w r , e + f"wr, e - 2 l K f'^J 

3 

We now evaluate / C\U(kz). Arguing as for the derivation of (23 1 we find the following identity 

/ U(kz) [fhU(kz) + 2fh'U{kz) + 2f'hk'VU{kz) ■ z] = 0, 

Jul*- 1 

which implies 

/ dU(kz) = fh 
On the other hand, using integration by parts, we obtain easily 

/ C 2 U(kz) = f C 3 U{kz) = 0. 

Turning to C4, we recall that 

C r w r = -2(f) 2 hU(kz)(H 7 $ + z) - h(X7 N V, $ + z)U(kz) + - ^ d jgil hkdjU{kz) 

3 

Therefore, after some computations we deduce 
/ C A U{kz) = -\( (£ r w r )m i -2fcVi? J [ d j U{kz)w i -{p-l)h p - 2 I t/(fe) p " 1 u; r 

Recalling the expressions of £ r and Ci and integrating by parts we then get 

/ C A P{kz) = -\! (jC i w i )w r -2kY'H j I djU(kz) 



(('; 
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From the definition of uii and using some cancelation we then find 

1 



(C 4 + C 5 )U(kz) 



2hf J Wl U(kz) 
2hf [ w r>e U(kz) 



2kJ2H j J d J U(kz)w l 

j R " 

-2kJ2 Hj J d j U(kz)w l . 



To evaluate the last integral we need to use the explicit expression of w r . e : recall that we have 

h[(v N v^} + 2(f) 2 (ii^) 



(p-l)hP 

Therefore, adding all the terms and using some scaling we obtain 

fh 



' Y U(kz) + -^-\7U(kz) ■ z 
1 2k 



Ri, e U(kz) 



k n ~ 



k n - 



l -d s Uh 2 f [(V N V, 4>) + 2(/') 2 (H, $> 



1 / 1 n-1 
k n+1 \p-l 4 



U z 



From the Euler equation (34 1 it follows that 



(V N V,E m ) + 2(f) 2 H m = ?-^hP- l H m . 



Therefore after some manipulation we find 



f7 2 8,(H,$). 



Contribution of Rj. e .. f, • M ultiplying i?i. e ,/ 1 by hll(kz) and integrating, recalling the expression of Wi 



determined in subsection 



3.2 



(wi.o = - Y,j%f'hz U(kz) and w ife = p -^- fh'\z\ 2 U(kz)), we obtain 



Ri, e ,h U = 2hh fi U z + 2h z f[k' \7U-zU 



i .„;ir>iJ i ^ .rr , P —f'2f[hh' I \z\ 2 U 2 



+ f"h 2 U 2 +4h 2 f 2 f{' UU + 4hf'f[ Ud s (hf'U) 



UU - 



+ 2h 2 f"f'f[ 



h p h'f' 2 f[ / U p \z\ 2 U 



n — 1 
2fc™ 



We will need the following observations 
/ VU(kz) ■ zU(kz)dz = — V / djU{z)zjU{z)dz 

which imply 

(56) / UUdz = [ t * _ 1 U 2 (kz)+^-\7U(kz)-zU(kz)dz 

jR»-r 7 R „-i (p- ljftP 2fc 



t/ 2 (z)dz 



n — 1 
2fc™ 



2(p-l)*»+ 1 V 1 
/ \z\ 2 \7U(kz) ■ zU p {kz)dz = — - V / z?d j U{z)z j U»{z)dz 

= - ^ E /> + i)^u^(z)dz - -L £ £ pjdMttewWz 



U 2 {z)dz- 



U 2 (z)dz; 



n+1 



\z\ 2 U p+1 (kz)dz - p / |z| 2 V?7(fcz) • zU p (kz)dz. 
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From these we get 



\z\ 2 VU(kz) ■ zU p (kz)dz = ~ U + \ [ \z\ 2 U p+1 (kz)dz = -- " ' 



\z\ 2 U p+1 (z)dz 



(p+l)kJ R J (p+l)k n + l 
Moreover one can easily see that 

2hh'f[ f U 2 dz + 2h 2 f[k' f VU ■ zUdz + f{h 2 f U 2 dz = djfh 2 f[ f U 2 (kz)) dz 

h 2 f[ 



and 



Ik 



Ahff[ Ud s (hfU)dz = 4^{h 2 f 2 f[ UU)dz + 4f'hWf[] 



U\z)dz 



UUdz 



±f' 2 h 2 f[ / UdsUdz 
ii"- 1 

4ds(h 2 f ,2 f[ [ Uu)dz + 4fhds[hf'f[] [ UUdz 



4f 2 h 2 f[k' 



(p-l)hP- 1 J K n-l 

V' 2 h 2 fX I \z\ 2 \VU{kz)\ 2 dz 



U{kz)VU{kz) ■ zdz 



On the other hand, by (64) (see Section 5.1l we have 

P ^ £/ l /" I.,l2rr2j_. P — 1 

P+l' 



JM.n-1 



f 2 f[hti / |z| 2 t/ p+1 (fcz)dz 



Hence with some easy manipulations and by ( 65 1 one finds 



V -^f 2 f[hti \ \z\ 2 U 2 dz !/ ' / r-/-i.i2 f -./. .)f^.-f 



|z| 2 |VJ7(fcz)| 2 dz 



( P -lKn-l) //2/ ^ 
(p-l)(n-l) 



(i 



« + 1 7m»-i 



|2| 2 |V/7(fcz)| 2 dz- 



p-1 
p+l 



\z\ 2 U p+1 (kz)dz 



4k r 



^f'Zfihti f U 2 (z)dz 



Collecting the above expressions, and letting Cq — L n _i U 2 (z)dz one proves that 



ffc e f, C/dz = dg I ^4 I - 29f f 2 h 2 f[-, — - — r 



k 

f'2 flu' 



(p-l)(n-l) 2 

4A ,n+l 



V rz fxhh'-, ^r-r + (n- l)f ,2 f[^ T + f'f'fih 2 -, ^--j 

1 (p-l)fc n+1 v ' J1 k n+1 J J1 (p-l)k n+1 



2 f i r 



1 - 2/' ; 



(p- l)fc 2 



using (23 1, we have that af' 2 h' — f'f" = 0, which implies 

'h 2 f[ r 



(p - l)fc™+ 



— hf[(af' 2 h'-f'f"); 



h 
Co 



R heJl U(kz) 



^ 2 /I 
(p - l)fc n + 



1 - 2/' 2 



(p- l)fc 2 



— [(p-l)^- 1 -2^ 2 / l 2CT ] 
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Recalling the discussion at the beginning of this subsection, we want to find f± such that 

/ R i , eJl U(kz)+ I R he U(kz)=0. 

This is equivalent to 

h 2 .n 



(57) Tfi := 8, I fa _^ ra+1 [(p - l)^ 1 - 2aAW] 

Hence we get 



(p-l)t 
2A(p- l)k n+1 



P-1 

- l)/^ 1 - 2<7.4 2 /l 2ff + 2 V 20 



1 (H,$) 



2^(^-l] ( V(H.«I»> 



(p-l)k n+1 
\ P -\)hP+ l -2oA 2 h 2 °+ 2 ' 



where c is a constant to be chosen so that L f[ds = 0. Noticing that, by (31 1, we have 



- \)hP^ - 2oA 2 h 2 °+ 2 OA 



the required condition becomes 

'p-1 
26» 



i )2A(p-iy- 



I Aah 



cr-l 1 



8A 



h a ds 



As one can easily check from ( 29 1 and ( 34 1 , c coincides with A' and therefore we have in conclusion 

(p- l)fc" +1 



(58) /( 



2.A(p- l)fc" +1 

(p - l)/jP+! - 2(7.4 2 /l 2 



P-1 

20 



1 (H,Q)+A' 



( P -l)hP+ 1 -2aA 2 h 2 °+ 2 ' 



4.2 Projection onto the kernel of C r 



Similarly to the previous subsection, we need to annihilate the projection of the 7?'s in (551 onto the 



kernel of C r . This corresponds to multiplying the error terms by d m ll(k-), m = 1, . . . , n — 1, integrating 
over M n_1 and taking the real part. As before, we are left to consider only two terms: R ro and R r o j 1 . 
We beg in b y multiplying R r o by d m U and integrating. Recalling the expression of Wi determined in 
{wi,o = ~£j ^'■fhz j U{kz) and w he = ^fh'\z\ 2 U(kz)), we obtain 



Subsection 



3.2 



where 



(Rr >0 + R Tt0jl )d m U(kz) = Ai + A 2 + A 3 + A A + A 5 + A*, 

-1 

A 1 = 2(f) 2 ^h [ z m U(kz)d m U(kz) + &^hk [ (d m U{kz)f 



a, = 



2ff"h J i z m U(kz)d m U(kz) + 2(f') 2 h' i z m U{kz)d m U{kz) 
+ ff'h [ z m U(kz)d m U{kz) + 2h'k [ (d m U(kz)) 2 + 2hk' [ {d m U(kz)) 2 
+ 2{f) 2k -h [ z m U(kz)d m U(kz) ; 

fc JR"- 1 



A* = $L 



2{f ) 2 hk' [ z m d m U(kz)VU(kz)-z + 2hk'ky2 I zid 2 ml U{kz)d m U{kz) 

|^| 2 (c/(^))^ m a m f/(^) 



+(/')-*■'/, / j^(a m f/(fo)) 2 + i(p-i)(/') 2 / l ^ 
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Here we used the fact that ^ = v -^~ \ ■ We have next 



A 4 = 2(f) 2 H m [ z m d m U(kz)w r . e + 2(f) 2 (H,<S>) f w r . d m U(kz) 

d m w rie d m U(kz) -p{jp- l)h p ~ 2 / U(kz) p ~ 2 w r . e w r . d m U(kz) 



H" 



+ / (V A V, W rte Z + W r:0 Q}d m U (kz)] 



A* = 



and 



-{f) 2 hJ2^idf mgil f z m U(kz)d m U(kz) + 8(f) 2 h(H,^)H m [ z m U(kz)d m U(kz) 
x ii"- 1 Ji"- 1 

hk(H,$)H m [ (d m U(kz)) 2 + hk 2 V df s g tj ^ f d 2 3 U{kz)z s d m U{kz) 
hk V df t g tm $i ! (d m U(kz)) 2 + hk(H, <S>)H m f (d m U(kz)) 2 

\hk[Y^d 2 migil <S>i\ f {d m U{kz)) 2 + Y,d 2 mj V$ j h [ z m U(kz)d m U(kz), 

Z \ l J JR™- 1 ■ JR"- 1 



A 6 = 2ff[[ w r , d m U(kz) + 4ff[hH m f Uz m 8 m U(kz) 

JR"-1 JE"- 1 

- 2p(p-l)f'f{h?- 1 [ U p ~ 2 Uw r , d m U(kz)+2f'f[h(V N V) m [ z m Ud m U(kz) 

JR™- 1 JR™- 1 

+ 2H m f'f[hf d m Ud m U(kz) + AH m {ffff[[ z m d m U(kz)hU. 



Integrating by parts and using the above relations we obtain 
h 



At = 



k" 



U 2 {z)dz^^l % ^- V —^-h° 
y ' m k n 29 



p-\ 



U 2 (z)dz. 



Similarly, for A2 after some integration by parts and some rescaling one finds 



A 2 = 



dz 



k 



U z (z)dz } . 



In particular, recalling the identity (45) one then gets 

&m fp-1 



A 2 = 



(hk)' - 



Now we turn to the term A3 . First of all we can write 

n-1 

^ J ^ zid; nl U (kz)d m U (kz) = - 

1 



U 2 (z)dz. 



V/ z l d 2 ml U{kz)d m U{kz)=- V —± f (d m U(kz)) 2 -V / Zl d 2 ml U(kz)d m U(kz). 



Hence, using a simple scaling and again formula (45) we obtain that 
V / z l d 2 ml U(kz)d m U(kz) = - {n ~ l){lp - l) 



46>fc r ' 



U 2 {z)dz. 
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Regarding the other terms in A3, we write them in radial coordinates, obtaining 
P\2 hU roo l{p-l){ffhP k' 



3^^u,„_ 2 r r n {U'{kr)fdr + 
Jo 



2 (n-1) fc 



/>OG 

„_ 2 / r n+1 U{kryU'{kr)dr 
Jo 



where w„_ 2 is the volume of S n 2 . Using a change of variables and integrating by parts (recalling the 
relation k 2 — we then find 



^- 2 (/') 2 k ' h [afrW* K» + i)(p-i) 



1 fc"+! 



2 p+1 



l [/ p+1 (j-)dr 



Using the formula in the appendix we find that this expression becomes 

n+l{f) 2 k'h 



2 k 



Yk'h f tt2 . . , 



Therefore we also obtain that 

"n + 1 (.f) 2 ^ 



2 



/ C/ 2 (z)& 



Finally, after some tedious but straightforward computation one also deduces 



A 2 + A 3 = $' 



p-1 h' t2 h 2a ti 
-A 2 a 



2 fc™- 2 



^p- 1 - 



1 2 ^ 2 ^, 



2rr 



p-1 



t/ 2 



It remains now to evaluate the contribution of A4 and of A 5 . Regarding A4 we recall that w r satisfies 



where 



F=(ffhU(kz) (j2 d m9n(z m + ^m)) +\Y, d ^ii hk djU(kz)~Y,^^ + ^ hU ( kz )- 

\ m / j m 171 

Differentiating this relation with respect to z m we obtain 

C r d m w r = d m F + p{p - l)h p - 1 kU{kz) lp -' 2 d m U{kz)w r . 
Multiplying by w r , integrating and using the self-adjointness of C r we find 

\pip - \)hv- 2 [ Uikzy- 2 d m Uikz)w 2 r --^rl F iFd m w r - w r d m F) . 



We have clearly 



Therefore it follows that 



= / d m iFw r ) = / iw r d m F + Fd m w r ) . 
\pip- l)h?- 2 [ Uikz)P- 2 d m Uikz) W 2 = / w r d m F. 

I J R n-l tlK J m n-1 



Hence A4 can be written as 



/ 



W r 







+ 


-f 


k Jr 


+ 


s 


+ 


h 2 k 2 



w r 



k J2 H l d 2 ml Uikz) + if'fd m Uikz) J2 Hl ( z > + 



(f'yU(kz)H m - kd m Uikz)(V N V, z + $) - hUikz)(V N V, E m ) 



R"- 1 
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After some cancelations and some integration by parts we find the following formula 



A 4 = -2k £ H l [ w r d 2 ml U(kz) + V —^k 2 f 
j Js.n-1 2k9 J R 



w r U(kz). 



Using the symmetries of the integrals we find 

-2kH r ' 



w r AU(kz) 
1 



p-1 



26 
p-1 



n - 1 



29 



n — 1 



n- 1 . 

1V re U{kz 

From the explicit expression of w r and some integration by parts we obtain 
1 ( 1 n-l \ f UP +i_ 



w r U{kz 
1 



U p (kz) 



1 



71 — 1 



n - 1 \p- 1 2(p + l) 



1 



n- 1 



p-1 



1 



P-1 

29 



Using (451 we find that / R „_ 1 U p+1 = 1 



(n-l)(p-l) 
29 



Ln-i an d therefore it follows that 



Next we turn to A;. First of all we consider the two terms 



TP 



B 1 = hk 2 V 9f s9ti $i / d 2 j U(kz)z s d m U(kz) + hky2df t9tm ^i [ (d m U(kz)) 2 . 
i, s ,t,j J * n - r Lt 1-1 

Looking at the first one, by symmetry reasons the summands do not vanish if, either s = m and t = j, if 
t — m and s — j or if j = m and s = t. In the first case, we see appearing the second derivative of g u , 
which vanishes by our choice of the geodesic coordinates. Therefore we are left with the terms 

B 2 = hk 2 V dlg^t j d 2 mj U(kz)z 3 d m U(kz) + hk 2 V df t g tm ^ [ d 2 tm U{kz)z t d m U{kz). 



Integrating by parts one easily finds 

B 2 = -\hkY j d? jgmj $t f (d m u(kz)) 2 - hk V df t9tm $i f (d m u(kz)) 2 



therefore it follows that 



B 1 = 0. 



We turn now to the other terms. Integrating by parts and using (451 one deduces 

-{ffh^idLm [ z m U(kz)d m U(kz)- l -hk[Y J d 2 ml g 1 i^i\ I (d m U(kz)f 



A 2 h 2a - h p 



-iP-1 
26 



Hence after some integration one finds 



20 



/ u 2 -^ya 2 v^ [ 



f/ 2 



4— (H,$)iT' 

k n 



-lP 1 _ A 2 h 2cr 
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Finally we turn to Aq: wc first evaluate the terms involving w r , , whose explicit expression is not 
known, but which can be handled via some integration by parts. Differentiating the equation C r U = —U 
with respect to z m we find that 

c r {d m u) = -d m u+p{ P -i)h p - 1 up- 2 ud m u. 

Therefore, integrating by parts and recalling the definition of uv, we obtain 

2/71 / w r . d m U{kz) - 2p(p - Vjf'Ah?- 1 f U p - 2 Uw r , d m U(kz) 

W ri0 £ r (d m U) = -2/'/i / dmUCrWrfl 



= -2/71 

= 2ff[hH 



dxU ( kdxU + V —^hP- l z x V 



Using ( 23 1 and ( 34 ) , we can also combine the last term in the second row of ^6 with the last one in the 
third row to obtain 

2ff[h(v N vy 



z 1 Ud 1 U{kz) + m m {fYff l h / z x dxU{kz)U 



= 2/71/J 



V —^W- 1 - 2A 2 h 2 ° 



z 1 Ud 1 U{kz) + 2A 2 h 2a H m 



z&UikzW 



2ff[h^—^h p - 1 H m 



z 1 UdiU(kz) 



With some manipulation, the sum of the last two formulas gives 

p- 1 h?- 1 



2ff[hH m 



Collecting all the terms in Aq then we get 



dxUdxU- 



2ff 1 hH m 2k 



dxUdxU- 



p - 1 hP- 1 



uu 



With some integration by parts, one finds that 

1 



d 1 Ud 1 U = 



n — 1 



U(-AU) = 



UU + 2 



n — 1 



UzxdxUikz] 



U(U P - U), 



and Ln-i Uz\d\U{kz) — — ^ J R „-i U 2 , and therefore the last quantity becomes 



2ff[hH m 



n — 1 



Using (|45| one also finds 

U{kz)U p {kz)dz 



p-1 



1 



UU ■ 



n — 1 



UU P - 1 



n — 1 



1 

¥\p-l 2(p+l 

i e p+i 

¥ (p+l)(p-l) 9 



U p+1 (kz)dz 

1 

k z p — 1 



U 2 (kz). 



From this formula, (56 1 and some manipulation we then get 

2a(p + 1) 



2ff[hH m 



2ff[—H 



2(n -!)(>- 1)6> (p- l)(n-l) 



U 2 (kz) 



p-1 
29 



1 



U 2 {z). 
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Now, recalling the expression of f[ given by (58 1, the total projection onto the kernel of C r is 



Co 



1 h 

2 fc™ 



2. 2cr 



h p - 



26> 



4— (H,$)iJ m 



fc™ 26* 



(p- i) 2 /i 



-iT n (H,$) 



p-1 



, P-1 V 



26> 2 /fc"- 2 

/if" 1 - ^rh 2 ° 



2A 2 o 



P-1 
29 



- 1 



h {2A) 2 {p-l)h' J+ -P- 1 
k"(p-l)h - 2aA 2 h< 7 



2 k n 



n^)H m - 



p-1 

a-9\ h 2AA'(p-l)h a+p - 1 



26 J k n (p- l)h e - 2oA 2 h? 



where Co — / R „_i U 2 . Now, recalling the definitions of a and 0, choosing fx as in (58 1, with some 
calculation we find 



p-1 1 



29 hk 



Co\-[h 6 - 



2 A 2 9 
~1 



(Rr t0 + R rt0 j 1 )d m U(k(es)z)dz 



ir } <i»:: ; - ,9 1 /," 1 - ) h'<$>' m + -^-h-°((v N ) 2 v)$ r , 



p-i 



p-i 



(59) 



2A 2 9_ 



p-1 yiZ^Vjmyn; I 1 [(p- - 2o\4 2 /i< T ] 



-(p - 1) (3 + f) h 29 - l^Mlh 2a + 2A 2 (5a + W)h e +° 



2A 2 oh° 



We notice that the operator between brackets coincides precisely with the one in ( 36 1 , corresponding to 
the second variation of the reduced functional which we determined in Subsection 



231 

Remark 4.1 According to the considerations in Subsection \2.4\ to every normal variation 0/7 it corre- 
sponds some variation in the phase due to both the variation of position and the variation of the constant 
A. Recall that the phase of the approximate solution is the following 



1 s 1 

F e = -f(es) = - 
e e 



fdl. 



Differentiating with respect to a variation v (see (12 )) we obtain 
d „ 



dv 



A v h° 



1 



Acrh 



1 fdh 
dA 



•At, 



dh dV\ 1 1 



dV dv 



2eJ 



Recalling formula ( 32 I we find 

dv £ e 2 A v 



i-i^h 1 
dA e 



Aah 



a _ x dhdV 11 
dV~dv + 2e J 



d„ gil Ah° 



d v9ll Ah°. 



Therefore, when we take a variation v 2 of"/, this also corresponds to a variation of the phase of j^F £ 
Notice that multiplying the horizontal part by hd m U corresponds to adding a variation of — ^r-. 
Hence integrating by parts we get 



A 



' P-1 
. 29 



- 1 



Mi 



P-1 a ,g-i9h , , „_t dh dV 1 
-2A A -- H dA- Aah dVdvr^2 Ahd ^ 
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4.3 Proof of Theorem 11.41 



The proof of Theorem 1.4 can be deduced with an iterative procedure, adding higher order corrections (at 
any arbitrary order) to the above approximate solutions. This method has been used for other singularly 
perturbed equations, and is described in detail for example in [35], [35], [32], [33]: hence, we will limit 
ourselves to a formal proof, since rigorous estimates can be derived as in the aforementioned papers. 
For meN we consider an approximate solution of the form 



1pn 



h(es)U(k(es)z) + e J (w r .j(es, z) + iwij(es, z)) 



y-®(es). 



Here y stands for normal Fermi coordinates as in Subsection 2.1 while we have set 



/(«) = / + Xy7i(«); 
J'=l 



where fj, <5>j are smooth real functions and normal sections respectively. We then write 



-Ag c tp m . e + V(ex)ip 



IP— 1 



] +o(e m ), 

v j=o 



where TZj are error terms depending on M, V, p, $ and /. 

In Lemma 3.2 we showed how to choose uv.i and Wi t i so that 1Z\ vanishes identically. In the previous 
two subsections instead we proved that IZ2 can also be canceled provided /1 satisfies ( 58 1 and <l>o (we are 
using the above notation) satisfies Z^o = 0. By the invertibility of 3 it is indeed sufficient to take both 
/1 and <&o identically equal to zero. 

Turning to higher order terms, we will find that the coefficient of e 3 will be of the form (up to the 
phase factor) 

C r w 3tr + i£iW 3it + G 3 (es, z), 

where Gsies, z) is an expression depending on V, u>i jr , wi^, W2, r ,W2,i, f, H and $1. As before, we will 
find that the above expression can be made vanish provided ji and $1 satisfy are periodic solutions of a 
system of the form (see (57 1 and ([59]) ) 



Tf 2 



( P 26 



.3.2 



l)^H,*i) + W 3j i in [(),£]; 

in [0,ij. 



Here W3,i and Wa^ are smooth functions of s independent of f% and $1. Again by our assumptions on 
7 the above system can be solved in fa, $1, and solvability up to order e 3 can be guaranteed. For higher 
powers of e one can proceed similarly. 



Remark 4.2 We stated Theorem \l-4\ for general powers in e for expository reasons. Indeed, for our 
purposes in l^OJj we will need approximate solutions up to order e 3 only. However, an accurate analysis 
of the error terms will be needed. 



Remark 4.3 If we multiply the operators 3 and T (see (36 1 and ( |57[ ) ) by h(s)k(s) and h(s) respec- 
tively, they become self-adjoint. This fact will be used crucially in the second part I4(rf , see in particular 
Subsection 2.3 there. 



5 Appendix 



In this appendix we collect some technical results: some integral identities first, and then the derivation 



of the second variation formula ( 36 1 . 



36 



5.1 Some generalized Pohozaev identities 

In this Subsection we derive some identities which are useful in the above computations. They follow 
from a simple integration by parts, but the proof is rather involved so we give some explanation here. We 
recall that the function U in (21 1 is radial, and so it satisfies the following ordinary differential equation 

(60) 



u" _ !L_V + u = u p 



Differentiating this relation with respect to r we obtain 
(61) - U 1 



— u" + + u' = P up- l u' 



We multiply (60 1 by r U" and integrate by parts to get (recall that we are in R n ) 



U'r n+2 U"'dr + (n + 2) J U'U"r n+1 dr + (n - 2) / U'r n+1 U"dr 
+ (n-2)(n+l) / {U') 2 r n dr 



(U'Yr n+2 dr-(n + 2) / UU'r n+1 dr 



p (U'yU p - L r n+ 'dr-(n + 2) U p U'r n+L dr 



Using (61 1 the above identity simplifies as follows 

/>oo 

(62) (n-5) / r n (U') 2 dr + (n + 1) / r"C/ 2 dr = 2 



n+ 1 



r n U p+1 dr. 



Similarly, if we multiply (60 1 by r t/' and integrate again by parts we find 



Ur n+L U"'dr + (n + 1) / UU"r n dr + (n - 2) / Ur n U" dr + n(n - 2) / UU'r n ~ L dr 
Jo Jo Jo 

UU'r n+1 dr-{n+l) U 2 r n dr = -p UU'U p ~ 1 r n+1 dr — (n + 1) / U p+1 r n dr. 



Using (61 1 the last identity simplifies as 
(63) 



(n-1) 



(Uydr= / r"(C/')^r + 



l U p+1 dr. 



From (|62| it then follows that 
(64) / r n U 2 dr 



P+Uo 



(" - 5) 
n+1 J 



r n (U') 2 dr. 



If we plug this identity into ( |63| we get 
(65) 



n-1 



6 



n+1 



\z\ 2 \VU\ 2 dz- 



P-1 
p+ 1 



|0| 2 [/ p+1 dz. 



5.2 Second variation 



The aim of this Subsection is to prove formula ( 35 1 for the second variation of the length functional stated 



in the beginning of Subsection 2.5 



37 



Proof of (351. Differentiating (28) twice with respect to t\,t% ( a * (ti,t 2 ) = (0, 0)), recalling our notations 
in Subsection 2.5 taking into account (19 1 and (20 1 we obtain 



(66) 
where 



Si + £ 2 = 0, 



A! 2 ah a - 1 ^(V N V,V)ds+ I Aa{a-l)h a - 2 



Aah° 



dV 
d 2 h 
dV 2 



dh 
dV 



(V"V,V}ds 



(V iX V, W) + 



d 2 h 

dVdA 



A 



(V iV V",V)ds + / Aah 



dV 



(iy N ) 2 v) [v,w]ds 



A' 2 J h a (H, V)ds - J Aah*- 1 ^ (W N V, V) (H, W)ds 

(H,V)ds; 



Aa \ h 
'7 



Ah" 



AA[a / h 

J -v 

A[a 



j j,m 

OA 

d 2 h 



V 3 W 71 



dh 

ds + A^A'zo- j ft'-^dS 



AA'[ 2 a J h^^ds + AA'Ma-l) J h° 



7 

cr-l 

7 

ft-" 1 



dh 
dA 



ds 



dAdV 



(V A V,W) + 



■4, 



Then the second variation of the energy is given by the following formula at (ti,t 2 ) = (0, 0) 



dhdt 2 



6(6 — l)h 6 



dh 
dV 



(V"V,V)ds 



6h° 



6h° 



_i dh 
dV 
d 2 h 



6h tt - L ^ ((V N ) 2 V) [V,W]ds 



dV 2 



(v Jv y,w) + 



d 2 h 



A 



dV K ' 1 dA 2 



dVdA" ' 2 

(H,V)ds 



<V"V,V}da 



a; / 0(0 - i)h 

Jo 



9/i 



-2 



r/.s 



9/i 



9/1 
&4 



r/.s 



d 2 h , M , S 2 /i ,, 



c/.s 



OAl^-'^ds- 
12 dA 



9h e - l ^{V N V,V) + dAh e 



dA 



(H,W). 



Now some cancelation will occur. We plug in the value of A'{ 2 from (66 1 into the last equality to obtain 

d 2 E E (u lPtit2 , Atit2 ) 



dt x dt 2 



= £ l + £ 2 + E 3 + E A + £ 5 + £ 6 + £ 7 + £ & + Ec>, 
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where the terms are given by 



£1= It 
Jo 



P- 1 io 



^V'lV' ^/,';,;„.V-'1V 



£2 = 



+ 



= 0(6 l)h e - 2 §^ [(V N V, V)A' 2 + (V N V, W)A[] 
+ l L 6h^J^[A' 2 (V N V,V)+A[(V N V,W}] 



^[^<V"V,V>+4<V W V;W>] 
2 ^j\a(a-l)h^^ [A[(V N V,W)+A' 2 (V N V, 



P 

P^iJo^W WV N V,W) + A> 2 {V N V,V)] ; 



V)] 



ds; 



t ^ S [A>1 (H ' W> + A ' 2 (H ' V>] + p~l t h<7 [A>1 (H ' W) + A ' 2 <H ' V>] 

AO f f)h 

\ A J ^ <T_1 ^[^i< H .W> + ^2< H .V>]; 



&4 



•4^0(0-1) / h e 
Jo 

\ A ^l h dA-—l AA ^- l) L h [OA) A 



L rP-h 

OA' A' h - 1 

yA,A 2 h dA2 

L 



p 

^AA> o t h°-^A> ■ 
P-1 1 Jo dA^ 2 ' 

^5 = ^%^- 1 ^((V^) 2 y)[V,W]-^^^/ l - 1 |^((V w )V)[V,W]; 



Jo 



eh 



0-1 



fd 2 h 

OV 2 



(V^V.VXV^V.W) 



2A " fL Aah^^(V N V,V)(V N V,W); 



2A6 r 



^ = -J q L ^ W) (H, V) + ^Aa £ h*- 1 (V N V, W) (H, V) ; 

£ 8 = - £ °h°-^ (V N V, V) (H, W) + ^Aa £ h"- 1 gj (V N V, V) (H, W) ; 



£9 = 



^-l)^- 2 -^a(<7-l)/i CT - 2 
p-1 



dh_ 
dV 



(v iV t/,v)(v iV y,w> 
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Now we will simplify each of these terms. We get immediately 

-l r 



Si = 



2A 2 8 
p-1 



Recall the identity ( 32 ) 



he _, dh _ 2A 



l dh' 



ds. 



dA P -l\ h ' + Arh ' OA 



differentiating with respect to V and A we get the two formulas 

dh 
W 



m e-a dh dh | h o-i d2fl _ 
' dV dA dAdV p-1 



2A ( i(j i dh A . ^ ua _ 2 dhdh A j u \ d 2 h \ 

oh" 1 — + Aa(a - l)h a 2 — — + Aah a 1 ] 



dVdA 



dAdV J 



1 3ah a 1 — + Aa{a - l)h a 2 



p-1 
2 

/ 

p-1 



dA 



dA 



Aah 



dA 2 



From these it follows immediately that 



20 f dh 
82 = 83 = 0; £4 = AA r I Aah' 7 ' 1 — + h c 



'p-1 



dA 



By means of (31), the above identity (32 1 can also be written as 



//" '2.4/r' — 1 '"'' ,T j '> l, ;iT 



dv P -i\ A ° h °-^W + h ° 



from which it follows that 
(67) 

This formula implies 



e _ x dh_ _ 2A 2 a t dh _ 1 

dv p-i dv p-i • 



s 5 = — - / ((v"yv)[v,w]h- 

p-lJo 



£7 = -j^jj j\v N V,W)(H,V)h-°; 
Therefore in conclusion we get at (ii,t2) = (0,0) 
9 2 E £ {u^ tit2 . Atit2 ) 



dhdt 2 



h° - 

p-1 



(p- 



T) j\v N v,v){n,w)h-°. 



3,m 



ds 



_! 2A*e_ K ,_ 1 \ d 2 h 



eh»- L - 
\ p-i 

L 



(v A v,v)<v A v;>v) 



p^i i> {(( vAr ) 2F )[ v ' w ] - (v w v,w)<h,v) - (v iV v;v)<H,w)j/r CT 



1)^-2 _ M!^ (cr _ 1)h a-2 
p- 1 



29 f ( UCT _ X dh ur 
P-1 Jo V 9 A 



dh 
dV 



(v JY y,v)(v JY i/,w) 
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Differentiating (67 1 with respect to V, we obtain 



[(p - I)//" 1 - 2aA 2 h°- 1 ] + [(p - 1)(0 - I)//" 2 - 2a(a - l)^ 2 ^ 2 ] 



.1 9ft, 

9V' 



which yields 



2A 2 

( 

p-1 



Oh 9 - 1 - ^^ah*- 1 



2 

dV 2 

L 



,yh (v N v,v)(v N v,w) = - Vn 



dh 



P-1 Jo dV 



(y"v,v)<y"v,w)h 



I) 



2a{a-l)A 2 e 2 
p-1 



l)h' 



8-2 



dh 
dV 



(V JV U,V)(V iV U,W) 



Collecting the above computation together with some further cancelation, one finds, at (ti,t 2 ) — (0,0) 



dhdt 2 



p-1 



ds 



P-1 Jo 
ad 



{((V N ) 2 V)[V,W] - (V N V,V)(H,W) - (V N V 7 W)(U,V)}h^ds 



P- 1 J a 



dV 



29 
P-1 Jo 



This proves formula ( 35 1 . 
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